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TAP 1: A FINITE ELEMENT PROGRAM FOR STEADY-STATE THERMAL 

ANALYSIS OF CONVECTIVELY COOLED STRUCTURES 


By Earl A. Thornton^ 

SUI^MARY 

A finite element computer program (TAP 1) for steady-state thermal 
analysis of convectively cooled structures is presented. The program has 
a finite element library of six elements; two conduction/ convection 
elements to model heat transfer in a solid, two convection elements to model 
heat transfer in a fluid, and two integrated conduction/ convection elements 
to represent combined heat transfer in tubular and plate/fin fluid 
passages. Nonlinear thermal analysis due to temperature dependent thermal 
parameters is performed using the Newton-Raphson iteration method. Program 
output includes nodal temperatures and element heat fluxes. Pressure drops 
in fluid passages may be computed as an option. A companion plotting program 
(TAPPLT3 for displaying the finite element model and predicted temperature 
distributions is presented. User instructions and sample problems are 
presented in appendixes. 


^ Associate Professor of Engineering, Old Dominion University, Norfolk, 
Virginia 23508. 



INTRODUCTION 


TAP 1 (Thermal Analysis Program) was written in the course of 
research aimed toward the development of finite element methodology 
for the thermal analysis of convectively cooled structures. The finite 
element methodology and applications to several convectively cooled 
structures are presented in reference 1. 

The main body of this report presents: (1) the salient features of the 

finite element theory, (2) the computer program organization, (3) the finite 
element library, (4) the nonlinear solution algorithm, and (5) a companion 
computer graphics program TAPPLT. General directions for use of the programs 
are presented in Appendices A-C, Program input data and output are illus- 
trated with sample problems in Appendix D. 

FINITE ELEMENT FORMULATION 


Thermal analysis of convectively cooled structures includes coupled 
conduction and convective heat transfer in a region consisting of a solid 
structure and a moving fluid. The problem may be mathematically formulated 
in terms of the energy equations of the solid and fluid assuming incompressibl 
flow (ref. 2) . For steady-state heat transfer, neglecting viscous energy 
dissipation in the fluid, the temperature T (x, y, z) satisfies: 

->■ -> 

V • (Kg * 7T) = 0 (1) 


for the solid region and 


pCp V • VT - V . (|^ • 7T) = 0 


( 2 ) 


for the fluid region where and Kp denote the 
the solid and fluid, respectively; p is the fluid 
is the fluid specific heat. The thermal properties 
fluid (^ and C^) are temperature dependent. The 
(eq. 2) specifies the fluid motion as a function of 


conductivity tensors of 
mass density; and 
of the solid CK_) and 
velocity vector V 
the spatial coordinates 



and, in general, is unknown. Equations (1) and (2) must be solved subject 
to specified boundary conditions on the external surfaces of the region 
and appropriate continuity conditions at the solid/fluid interface. 

In the most general approach to thermal analysis of conv actively 
cooled structures the steady-state velocity distribution of the fluid must 
first be determined by solving the momentum and continuity equations of 
fluid flow. With the fluid velocity distribution known, equations Cl) 

(2) may then be solved simultaneously for the temperature distribution in 
the solid/ fluid region. 

In. TAP 1 a simplified finite element solution procedure is employed 
using a number of assumptions customarily used in practical heat transfer 
analysis (ref. 1) . The thrust of the assumptions is the elimination of 
the computation of the fluid velocity distribution. 

The solid region of the convectively cooled structure is represented 
by standard conduction/convection elements. Two conduction/convection 
elements, a rod and a quadrilateral, are available in the program. The 
fluid region of a convectively cooled structure is modeled by elements 
which represent convective heat transfer in the coolant passages. Basic 
convective finite elements were developed for the fluid to represent each 
of the terms in equation (2). The first term is represented by a mass 
transport convection element, and the second term is represented by 
surface convection elements with unknown fluid temperatures. Special 
integrated conduction/convection finite elements were also developed; (1) a 
tube/fluid element, and (2) a plate-fin/ fluid element. The basic convection 
elements and the integrated conduction/convection elements may be combined 
with the standard conduction elements for analysis of a variety of convectively 
cooled structures. Any of the elements may also be used independently. 

PROGRAM ORGANIZATION 


The organization of T.'\P 1 is based on the SAP family of structural 
analysis programs (refs. 5, 4). A flow chart of the TAP 1 main program is 
presented in figure 1. The main program consists of eight subroutines which 
are sequentially called in a normal program execution. These subroutines 
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Figure 1, Flow chart o£ program TAP 1. 
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process input data, generate plot files, assemble and solve the equations, 
print nodal temperatures, and perform heat flux calculations. Dynamic 
storage allocation is used to store all input data and large arrays in a 
blank common designated in the main program as A. At execution, the amount 
of blank common storage available is calculated in the main program from 
the JOB card field length. The amount of blank common is the only restriction 
on the amount of input data, i,e., there are no other limitations on number 
of nodes, elements and thermal parameters. 

Nodes (INPUT) 

The thermal system is described by a set of nodal points with unknown 
temperatures. A nodal point is described by a data card containing the 
node number, a boundary condition code (zero or one), the nodal coordinates, 
a generation parameter, and a specified nodal temperature if required. 

Nodal points which may be entered in rectangular Cartesian coordinates 
Y) z) or cylindrical coordinates (R, Y, 0). Input data for regular 
nodal patterns may be reduced by utilization of a nodal generation 
capability based on linear interpolation. All nodal point data are 
retained in core during the assembly of the element conductance matrices. 

Nodal point data are also saved on an auxiliary storage file if plots are 
requested. 

A boundary condition code of zero indicates an unknown nodal temperature. 
A boundary condition code of one indicates a specified nodal temperature 
which will be held constant during the solution. Heat loads and convective 
boundary conditions are specified as a part of the element input data. 

Elements (ELTYPE) 

Elements are entered into the program in groups which consist of a 
number of sequentially numbered elements of the same t>pie. There may be 
more than one group of the same element type. Data generation schemes are 
provided for all elements to reduce input data for regular finite element 
meshes . 

The input data for all elements follows the same general scheme’ 

(1) a control card for each element group, (2) a set of thermal parameter 



cards, and (3) a set o£ element cards. For a linear analysis thermal 
parameters are entered as constants; for a nonlinear analysis table 
numbers are entered. Each element may have different thermal parameters. 

Element conductance matrices, heat load vectors and heat flux matrices 
are computed as the element data cards are read. These matrices are stored 
sequentially on files for later use in assembly of the system equations and 
in heat flux computations. For elements with more than one thermal parameter, 
the element conductance matrices are resolved into components, one for each 
thermal parameter. For a linear analysis the conductance matrices are formed 
for the thermal parameters entered; for a nonlinear analysis the conductance 
matrices are formed initially for unit thermal parameters. Element 
connectivity data are saved on an auxiliary storage file if plots are 
requested. 

As element data are processed the system bandwidth is computed. 

Bandwidth is defined herein as the maximum difference between two connected 
node numbers plus one to account for the diagonal. The bandwidth is used 
later in the program to determine storage requirements for the system 
conductance matrix. For optimum program storage requirements and execution 
times the bandwidth should be a minimum. Bandwidth is determined by the 
user's nodal numbering scheme (see ref. 5). 

After all element data have been processed, the nodal coordinates are 
no longer needed, and the corresponding core storage area is used for other 
variables. The nodal boundary conditions are, however, retained in core 
since they are required later in the solution process. 

~ Thermal Parameter Tables (TABLES) 

For a nonlinear analysis thermal parameter tables are required. The 
input data consists of a control card for each table and a set of cards 
containing temperature and thermal parameter data points . Tne thermal 
parameter data are retained in core during the balance of the solution 
process. Ordinarily, the amount of core storage required for the tables 
IS small in comparison to storage required for the system matrices . In the 
solution process, linear interpolation and extrapolation are used in looking 
up thermal parameters. 
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Assembly of System Matrices (ADDSTF) 

The system conductance matrix and heat load vector are formed in 
blocks as shown in figure 2. Because of mass transport convection and the 
Newton-Raphson iteration process, the system conductance matrix may, in 
general, be asymmetrical. Hence, advantage of matrix symmetry could not 
be taken as in structural analysis. The number of equations in a block 
depends on the amount of blank common storage available and is computed 
in the program from the equation: 

NEQB = (MTOT - 2*NUMNP - 2*NUMTB - KENT) / C4’'MBAND3 (3) 


where 


NEQB - number of equations in a block; 

MTOT - amount of blank common storage; 

NUMNP - number of nodal points; 

NUMTB - number of tables; 

KENT - number of table entries; 

MBAND - bandwidth of conductance matrix. 

With the number of equations per block known, the thermal conductance matrix 
and heat load vector are assembled two blocks at a time by direct addition 
of the element matrices. During a nonlinear analysis the Newton-Raphson 
correction to the element matrices is performed as part of the assembly 
operation. In the assembly operation it is necessary to pass through the 
element matrices which are stored on a file. In order to minimize file 
reading (ref. 4), element matrices which pertain to the next several blocks 
are written on another file. This method, for large problems, significantly 
reduces the amount of file reading in the formation of the equation blocks . 
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Figure 2. Program storage of system matrices. 
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In a linear analysis the equation blocks are assembled only once; in a 
nonlinear analysis the equation blocks are reassembled for each iteration 
(see fig. 13 . 


Boundary Conditions (TEMPBC3 

In finite element thermal analysis with TAP 1 the only nodal boundary 
conditions required are specified temperatures (i.e., temperature 
gradients can not be specified) . Specified nodal temperature data are 
entered into the program with the nodal input data. Heat fluxes and 
convective boundary conditions are entered with the element data and are 
incorporated by the program into the system heat load vector. For a 
boundary with zero heat flux, no boundary condition need to specified; the 
heat load terras corresponding to zero heat fluxes are automatically taken 
as zero. 

The program handles the temperature boundary conditions using the method 
described in reference 6. Basically, this method consists of modifying 
the conductance matrix and heat load vector such that the size of the 
matrices is unchanged. The advantage of this approach is the ease of 
indexing the equations, i.e., the node numbers and equation numbers are 
the same. A disadvantage is that extra equations are carried in the 
solution process. For TAP 1 thermal analysis temperature is the only 
degree of freedom per node, hence the penality is not very large since usually 
only a small percentage of the equations have specified temperatures. 

Solution for Temperatures (TAPSOL) 

The general, banded, simultaneous equations are solved by Gauss 
eimination. The subroutine is based upon the banded out-of-core symmetric 
equation solver, BANSOL, presented in reference 7. Operations with zero 
coefficients are skipped. Matrix data is transferred into core two blocks 
at a time. If a sufficient amount of blank common is available to store the 
equations in two blocks or less, an in-core solution is performed. The 
basic limitation on the equation sol't^er for a given field length is that 
the number of equations in a block must be greater than the bandwidth. 
Normally, this restriction poses no problem and may be overcome by 
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increasing the field length or renumbering the nodes to reduce the 
bandwidth. Block size is automatically determined at execution. 

For some assemblies (e.g. , in series) of mass transport convection 
elements it is possible to obtain zero coefficients on the diagonal of the 
conductance matrix. Dependent on the boundary conditions, such a problem 
may cause the equation solver to stop with an error message to avoid a zero 
divisor in the Gauss elimination process. This difficulty can normally 
be overcome by renumbering the nodes so that a zero diagonal coefficient 
is filled in during the elimination process, A. zero diagonal coefficient 
will not arise in the integrated thermal/fluid elements for a nonzero 
convection coefficient. 


Heat Flux Calculations (FLUX) 

After the nodal temperatures are computed, element heat fluxes are 
calculated using element matrices previously stored on a file. Typical 
element fluxes calculated include, e.g., for the quadrilateral conduction/ 
convection element, conduction heat flux components at the element centroid 
and convection heat fluxes on the top and bottom surfaces and four edges. 

In' general, conduction heat flux components are positive in directions of 
the local element axes , and surface convection fluxes are positive into a 
surface. 

For the integrated thermal/fluid elements pressure drops are computed 
as a user option in the heat flux computations. Pressure drop computations 
include flow-friction and flow-acceleration effects (see ref. 8) . Pressure 
drops are computed for three user options; (1) constant density, (2) variabl 
density using a density -temperature table, or C3) an ideal gas. 

THE ELEMENT LIBRARY 

The library consists of six elements (fig. 5) for either linear or 
nonlinear theimal analysis. 



(a) Conduction/Convection 
Rod Element. 


(b) Conduction/Convection 
Quadrilateral Element, 



(c) Mass Transport Convection Element. 




(d) Surface Convection Elements . 



Figure 3. Element library of TAP 1, 
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Conduction/Convection Rod Element 

The rod element is based on the same assumptions as the NASTRAN rod 
element (ref. 9). A linear temperature variation is assumed between nodes. 
The element permits heat loading due to internal heat generation, prescribed 
surface heat flux or surface convection. The convection heat transfer 
coefficient and fluid medium temperature may be different at each node. 

Conduction/ Convection Quadrilateral Element 

The quadrilateral element is based upon an isoparametric formulation 
similar to the approach described for structural elements in reference 5. 

The term isoparametric means the same interpolation functions define the 
element shape and the element temperature distribution. The temperature 
within the element is given by 

4 

T (?, n) = E N. T. (4) 

where are the interpolation functions , 

Ni = j (1 - (1 - n) 

N 2 = I* (1 + €) (1 - n) 

CS] 

N3 = (1 + C) (1 + n) 

N4 = i (1 - 5) (1 - n) 

and T^ are the nodal temperatures . The quantities 5 , n denote the 
isoparametric coordinates for a unit square. The conductance matrices are 
computed for the element by integration in the ? , n plane; in TAP 1 
the conductance integrals are evaluated by the four-point Gauss quadrature 
rule of numerical integration (ref. 5) . 
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For rectangular elements, the conduction heat flux component varies 

linearly with y , but it is independent of x ; similarly the component 
varies linearly with x , but it is independent of y . Conduction 
heat flux components are always calculated at the element centroid. 

The quadrilateral element permits a laminated composite material. 

Each lamina 'is assumed to be orthotropic; input data for a lamina consist 
of a conductivity tensor, a material axis angle and the lamina thickness. 

An arbitrary number of lamina are permitted. For a nonlinear analysis the 
lamina properties are assumed to have the same temperature variation, i‘. e., 
an element is characterized by a single conductivity-temperature table. 

The element permits internal heat generation, prescribed edge or 
surface heating, and convection heat transfer on all four edges and 
the top and bottom surfaces. Convection coefficients and fluid medium 
temperatures may be different at each node. 


Mass Transport Element 


The mass transport element represents energy transport downstream due 
to fluid flow. The element represents the first term in equation (2) and 
is based on the following assumptions (ref. 1) : (1) the thermal energy state 
of the fluid is characterized by the fluid bulk temperature which varies 
only in the flo\if direction, and (2) the fluid velocity is represented by a 
mean velocity V which varies only in the flow direction. The basic input 
data IS the mass flow rate, m , and the fluid specific heat c 


The mass transport element has an indefinite, asymmetric conductance 
matrix [see ref. 1). As previously discussed (see Solution for Temperatures), 
some assemblies of mass transport elements may create zero diagonal terms in 
the system conductance matrix. 


The computation of the element conductance matrix does not depend on 
the coordinates of the fluid nodes. Thus, fluid nodal coordinates are 
arbitrary and are used only in plots. 
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Surface Convection Elements 

Surface convection elements fa quadrilateral and a triangle) represent 
the energy transfer between a coolant passage surface and the fluid. The 
heat transfer is based upon a convection coefficient for the fluid and a 
surface area of the passage. The surface area is computed as the product 
of the distance between wall nodes and an area factor supplied as input data. 
The fluid nodal coordinates are arbitrary and are used only in plots. 

Tube/Fluid Integrated Element 

The tube/fluid element consists of fluid within a thin tube of constant 
thickness and constant, arbitrary cross-section. The element has two 
fluid nodes I and J and two tube nodes L and K . The fluid node 
locations are arbitrary at a given flow section and are used only in plots , 

The following heat transfer modes are represented in the element; 

1. Axial conduction in the tube between nodes L and K ; 

2. Convection between the internal tube surface and the enclosed 
fluid (nodes L, K, and nodes I, J) ; 

3. Mass transport convection due to fluid flow from I to J ; and 

4. Heat transfer between the external tube surface and a surrounding 
medium which is represented by specifying a heat flux or the medium temperature 
and convective film coefficient. 

The convection area between the internal tube surface and the enclosed 
fluid is computed as the product of the distance between tube nodes and the 
input tube perimeter. The external heating is assumed uniform around the 
perimeter of the tube. The surface area for external heat transfer is 
assumed equal to the internal convection area. The temperature and convection 
coefficient of the surrounding medium may be different at each tube node. 

As user options. (1) the fluid convection coefficient may be modified 
for large temperature differences between the fluid and tube surface (ref. 8), 
and (2) fluid pressure drops may be calculated (see Heat Flux Calculations). 
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Plate-Fin/Fluid Integrated Element 

The plate-£in/fluid element consists of two walls (plates) connected 
by an internal fin. For convenience a single plain fin is shown in figure 
3; in practice other fin configurations (e.g.,, pin or offset fins) may be 
represented by using an equivalent thickness and surface area for the single 
fin. Fluid flows along both sides of the fins through an arbitrary flow 
cross section (shown trapezoidal for convenience), which may. vary linearly 
along the element. The element has 6 nodes: two nodes to represent the 

fluid bulk temperatures -(nodes N and K) and four wall/ fin nodes 
(I, J, L, and M) . The fluid node locations are arbitrary at a given flow 
section and are used only in plots. The following heat transfer modes are 
represented in the element: 

1. Two-dimensional conduction in the fin between the nodes I, J , 

L , and M ; 

2. Convection between the fin surfaces (nodes I, J, L, and M) and 
the fluid (nodes N and K) ; 

3. Convection between the wall surfaces (top nodes M and L; bottom 
nodes I and J) and the fluid (nodes N and K) ; and 

4. Mass transport convection due to fluid flow from N to K . 

The fin is modeled as an isoparametric quadrilateral element with 
surface convection to a fluid with unknown temperatures . Input data describing 
the fin includes its effective thickness and an area factor for convection. 
These quantities may be adjusted as input to permit the plain fin to 
represent other fin configurations. 

Convection between the wall surfaces and fluid is based on areas 
computed using input wall widths, the fin thickness, and internally 
computed distances between wall nodes. The flow cross-sectional area may 
vary due to a difference in passage height at the element entrance (I to M) 
and exit (J to L) . 

User options are available to. (1) modify the convective heat 
transfer coefficient for large temperature differences between the fluid 
and wall surfaces (ref. 8), (2) modify the fin convective heat transfer 
by an efficiency factor n which accounts for deviations in the fin 
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temperature variation from the assumed linear profile, and (3) compute 
fluid pressure drops (see Heat Flux Calculations). 

' NONLINEAR ALGORITHM 


The finite element fonnulation employed in TAP 1 leads to a set of 
nonlinear algebraic equations of the form 

[K (T)] {T> = {Q} (6) 

where [K (T)] denotes the temperature dependent system conductance matrix, 
{T} denotes the unknown nodal temperature vector, and {Q} is the system 
nodal head load vector. If thermal properties are not a function of temper- 
ature, equation (6) reduces to a linear set of equations which may be solved 
directly. If the thermal parameters are a function of temperature, the 
Newton-Raphson iteration algorithm is used: 

C7) 




C8) 


where denotes the system Jacobian matrix, and represents 

nodal residual heat loads . 

A key asstunption employed in TAP 1 is that thermal parameters are 
constant within an element. This assumption permits the nonlinear algorithm 
to be based upon one initial computation of element conductance matrices 
for unit thermal parameters. If a particular element depends on more than 
one thermal parameter, the matrix is formed by summing component matrices, 
one for each thermal parameter, TP . Thus a typical conductance matrix 
IS expressed as 


[f^3 = E 

m 


m 


C9) 
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where the siunmation includes all thermal parameters, TP^ , affecting the 
element, and [5^]m denotes a typical unit conductance matrix. For a typical 
element with N nodes the average element temperature is computed from 


T 


a 


E T 

A=1 


(103 


and a thermal parameter is looked up in the table using linear interpolation. 

The Jacobian matrix and residual load vector for a typical element are 
computed from 


J. . 

11 


TP*K. 

11 


^ 1 d(TP) 

N dT^ 
a 




R. 

1 



_ 'j'p-v 


N 


E 

A=i 



(iiD 


( 12 ) 


The quantity dCTP)/dT represents the slope of the thermal parameter 
curve. The formulations for the Jacobian matrix and residual load vector 
have the following computational advantages: (13 the equations are 

valid for all element types (i.e., rods, quadrilaterals, thermal/fluid 
elements, etc. 3, (23 element matrices need to be computed only once, and 
(33 the Jacobian and residual load vector have common operations. 

The Jacobian matrices and residual load vectors are computed for each 
element and assembled into system matrices at every iteration. The 
computations are performed as a part of the assembly operations in the 
subroutine ADDSTF called by the main program (see fig, 13. 

TAP 1 automatically uses zero nodal temperatures to initiate the 
nonlinear solution process. The iteration process is terminated when- 
(13 a specified number of iterations has been performed, or (23 the 
largest change in nodal temperature (expressed as a percentage3 is less 
than a specified value. For typical applications convergence has been 
obtained in from one to three iterations (i.e., two to four analyses 3 
using a convergence criteria of 0,1 percent. 
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PLOTTING PROGRAM 


A companion program, TAPPLT, is used to plot the finite element model 
and calculated temperature distributions. The program is based on 
the oblique orthographic projection program described in reference 10. 

The program includes options for plots of finite element models 
annotated with grid point or element numbers. Another option allows 
boundaries of an isolated portion of the model to be specified by cutting 
planes to permit detailed inspection of the selected region. Also, 
exploded views can be generated which separate the elements in a finite 
element model to detect the absence or presence of elements. Temperature 
surfaces, i.e., T = f (x, y) , can be plotted superimposed on the nodes 
of the model, or temperatures can be represented as vectors extending from 
the nodes. 

The program is limited to plotting elements with a maximum of four 
nodes so that the six-noded plate-fin/fluid element is plotted as two 
quadrilateral elements. 


CONCLUDING REMARKS 

A finite element computer program (TAP 1) for steady-state thermal 
analysis of convectively cooled structures has been presented. The program 
has a finite element library of six elements : two conduction/ convection 

elements to model heat transfer in a solid, two convection elements to model 
heat transfer in a fluid, and two integrated conduction/ convection elements 
to represent combined heat transfer in tubular and plate/fin fluid passages. 
Nonlinear thermal analysis due to temperature dependent thermal parameters 
is performed using the Newton-Raphson iteration method. Program output 
includes nodal temperatures and element heat fluxes. Pressure drops in 
fluid passages may be computed as an option. 

A companion plotting program (TAPPLT) for displaying the finite 
element model and predicted temperature distributions has been presented. 
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APPENDIX A 


PROGRAM DETAILS 

Computer and System Requirements 

TAP 1 and TAPPLT were written using standard FORTRAN IV and were 
developed on the CDC computer system at LRC. TAP 1 is almost system 
independent] only two system subroutines (JPARAMS and SIGN] are called. 

TAPPLT is dependent on the LRC computer graphics software for VARIAN and 
CALCOMP plots j but can be converted to other graphics systems (see reference 
10 for a documented source listing) . 

Storage Allocation 

Dynamic storage allocation is used by both programs. In TAP 1 all 
large arrays are stored in blank common designated as A; in TAPPLT large 
arrays are stored in blank common designed ZZZ. TAP 1 computes the blank 
common available from the job card field length and attempts to process 
the input data and perform a solution. The program terminates execution 
with an error message stating the additional storage required if insufficient 
storage is available. TAPPLT prints the required blank storage required for 
each execution; if insufficient storage is available, a system error message 
will result. The field length for both programs is problem dependent; the 
user may wish to estimate the required field length for any new problem 
attempted. 

TAP 1 Field Length Requirements . An approximate formula for the 
required field length (in octal) on the Network Operation System (NOS) at 
LRC is 


FL3 = 77,4003 T Ns (A- 1 ) 

where N is the additional blank common required which depends on: ( 1 ) the 

amount of nodal input data plus the storage required to process the element 
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type with the maximum amount of material data, and (2) the storage constraint 
imposed by the equation solver. 

For data imput the additional blank common required (in decimal) may 
be computed from 

Nio = 5*NUMNP + iM*NMAT 

where NUMNP is the number of nodal 
materials to be used for an element, 
parameters which can be input for an 


Element M 

Rod 3 
Quadrilateral 5 
Surface Convection 3 
Mass Transport 3 
Tube/Fluid 16 
Plate-fin/Fluid 16 


The second term in equation (A-2) is computed for each element group and the 
maximum value is used. 

For solution of the equations the minimum storage is determined by 
the constraint that the number of equations in a block [see eq. (3)] be 
greater than or equal to the system bandwidth. This requirement will be 
met if 


(A-2) 

points and NMAT is the number of 
M depends on the number of material 
element : 


Nio = 4*MBAND**2 + 2*NUMNP + 2*NUMTB + NENT (A-3) 

where NUMTB is the number of tables and NENT is the number of table 
entries. 

For optimum execution times it is usually desirable to have an amount 
of storage available larger than the minimum computed above (see Solution 
for Temperatures) . This guideline normally will insure that sufficient 
storage is available to process the input data. It may be helpful to note 
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that on the NOS a job card field length of 125,0000 will make available 
a blank common area of ll,000xo which is an adequate amount of storage for 
most moderate size problems (say, less than 300 nodes) . 

TAPPLT Field Length Requirements. A conservative formula for the 

■■ ■■■' ' • • • ' V 

required field length on the NOS is 

FLa = 33,000g + C5^NUMNP)8 . 

Auxiliary Storage Files 

TAP 1 uses 10 auxiliary storage files in a normal execution. The 
auxiliary storage files are identified in the table below. 

TAP 1 Auxiliary Storage Files 


File Function 

1 Element flux and pressure drop calculation data 

2 Element conductance matrices 

3 Temporary storage (TAPSOL) 

4 System conductance and heat load matrices 

5 Input data 

6 Printer output 

7 Temporary storage (ADDSTF, TEMPBC, TAPSOL) 

8 Nodal data for plots 

9 Thermal fluid modification data 

10 Element data for plots 
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APPENDIX B 


INPUT DATA FOR TAP 1 
General Setup of Input Deck 

The general setup of a typical input data deck is shown schematically 
in figure 4. A deck requires four basic data groups and a fifth optional 
group of data as follows : 

(1) A single heading card containing any desired title information; 

(23 A single master control card containing control values specifying 
various program options; 

(33 A node input deck containing nodal coordinates, the boundary 
condition code, and specified nodal temperatures, 

(4) An element input deck containing element data organized by groups. 
Each group consists of the following sequence of cards: 

(a3 a control card containing control values and a heading to be 
printed for the element group, 

(b3 a set of material property cards, 

(c3 a set of element cards; and 

(53 For a nonlinear analysis, the thermal parameter tables organized 
as a set of cards for each table. 

(a3 a control card containing control values and a heading to 
be printed with the table, and 
(b3 the data points in the thermal parameter table. 

Several problems may be solved on one program execution by placing the problem 
data decks in sequence. Plots can be obtained for only the last problem in 
a sequence. 


Input Data Cards 

Data cards are described in detail in this section. Input data is 
expressed in standard FORTRAN I, F, or A formats. Integers must be right 
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justified. The F format may be used to read real numbers expressed in an 
E format; however, numbers in an E format must be right justified. 

Any consistent set of units may be used. In the input data 
instructions which follow sample units are given for illustrative purposes 
only. 
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I . HEADING CARD 

(12A6) 


notes columns 

variable 

entry 

(1) 1-72 

HEDC12) 

Enter the heading information to be printed 
with the output 

NOTES/ 




(1) Begin each new problem with a heading card. 
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II, MASTER CONTROL CARD (7I5,F10.0) 


notes 

columns 

variable 

entry 

CD 

1 - 

5 

NUMNP 

Total number of nodal points in the 
model 

C2) 

6 - 

10 

NELTYP 

Number of element groups 

(3) 

11 - 

15 

NUMTB 

Number of thermal parameter tables (for 
nonlinear analysis) 

(4) 

16 - 

20 

NANA 

Analysis type 
.EQ.O; Data check only 

.EQ.l; Linear 

.EQ.2; Nonlinear 

C5) 

21 - 

25 

NDIAG 

Flag for diagnostic printing 
.EQ.O; No diagnostic output 

•GT.O; Diagnostic output 

(6) 

26 - 

30 

NPLOT 

Plot control code 

.EQ.l; Thermal model and input temperature 

plots 

.EQ.2; ^ Thermal model and computed 

temperature plots 

(7) 

31 - 

35 

NITER 

Maximum number of Newton-Raphson iterations 
(default .EQ.6) 

(8) 

NOTES/ 

36 - 

45 

TOL 

Convergence tolerance (default .EQ.0.1%) 


(1) Nodes are labeled with integers ranging from "1" to the total 
niiraber of nodes in the system, "NUMNP." 

(2) For each different element type (ROD, QUAD, etc.) a new element 
group must be defined. Elements within groups are assigned 
integer labels ranging from "1" to the total number of elements 
in the group. Element groups are input in Section IV, below. 

Element numbering must begin with one (1) in each different 
group. It IS possible to use ’more than one group of the same 
element type. 

(3) For a nonlinear thermal analysis, parameters are defined by thermal 
parameter tables. Thus for a nonlinear analysis (NANA .EQ.2) the 
number of tables should be entered. For a linear analysis 

(NANA .EQ. 1), the number of tables may be left blank or entered 
as a zero. 


26 



II, MASTER CONTROL CARD (concluded) 


(4) For NANA .EQ.O the program reads all input data, generating nodes, 
and elements as requested, and generates element matrices. Plot 
files are created for checking input data. Exit is made before 
the system matrices are assembled and the solution is performed. 

For NANA ,EQ.l the^thermal parameters are constant, and a linear 
thermal analysis is performed. An unsymmetrical set of banded 
equations is solved using Gauss elimination. 

For NANA . EQ. 2 the thermal parameters vary with temperature and 
are entered in tabular form. The equations are solved by 
Newton-Raphson iteration. 

For both linear and nonlinear thermal analyses conduction and 
convection heat fluxes are automatically recovered. 

(5) Diagnostic output may be obtained using this integer. This output 
typically consists of all element matrices, the eissembled matrices, 
and intermediate steps in the solution process. This option should 
only be selected for very small problems since a large quantity 

of data will be printed. 

(6) Two major plot options are offered: (1) For NPLOT .EQ.l plots 

of the thermal model may be obtained. In addition, plots of the 
temperatures input on the nodal cards may be made as either 

vector or surface plots (see plot instructions). (2) For NPLOT .EQ.2 
plots of the thermal analysis model may be obtained, and plots of 
the final computed nodal temperatures may be made in either vector 
or surface forms, 

(7) The Newton-Raphson iterative solution process will terminate when 
the number of iterations reaches the value NITER. Nodal 
temperatures are printed at each iteration, and element heat 
fluxes are calculated after the final iteration. The largest 
percentage change in nodal temperature will be printed at each 
iteration. 

(8) Convergence will occur if the largest percentage change in nodal 
temperatures is found to be less than the convergence tolerence. 
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III. 

NODAL POINT 

DATA (2(A1, 

I4),25X,3F10.0,I5,F10.0) 

notes 

columns 

variable 

entry- 

(1) 

1 

CT 

symbol describing coordinate system for 
this node]- 

EQ. ; (blank) cartesian (X,Y,Z) 

EQ.C; cylindrical (R,Y,0) 

C2) 

2-5 

N 

Node number 

(3) 

6 

I PR 

Print code 

EQ. ; (blank) normal printing 
EQ.A; suppress second printing 

nodal data 

C4) 

7-10 

IDCN3 

Boundary condition code 

EQ.O; No temperature specified 

EQ.l; Temperature specified 

(5) 

36 - 45 

XCN) 

X (or R) coordinate 


46 - 55 

YCN) 

Y coordinate 


56 - 65 

Z(N) 

Z (or 0) coordinate (degrees) 

(6) 

66 - 70 

KN 

Node number increment 

C7) 

71 - 80 

T(N) 

Nodal temperature 


Y 
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III. NODAL POINT DATA (continued) 

NOTES/ 

(1) A special cylindrical coordinate system is allowed for the global 
description of nodal point locations. If a "C" is entered in 
card column one (1) , then the entries given in columns 36-65 are 
taken to be references to a global (R,Y,0) system rather than to 
the standard (X,Y,Z) system. The program converts cylindrical 
coordinate references to Cartesian coordinates using the formulae: 

X = R sin 6 

Y = Y 

Z = R cos 0 


Cylindrical coordinate input is merely a user convenience for 
locating nodes in the standard (X^Y^Z) system, and no other 
references to the cylindrical system are implied. 

(2) Nodal point data must be defined for all (NUMNP) nodes. Node 
data may be input directly (i.e., each node on its o\«i individual 
card), or the generation option may be used if applicable (see 
note 6, below). 

Admissible nodal point numbers range sequentially from *'l" to the 
total number of nodes "NUMNP." Illegal references are: 

N.LE.O or N.GT. NUMNP. NUMNP must be the last card input. 

(3) The IPR variable is used to suppress a second printing of the 
nodal data. This would be desirable in the event that all nodal 
data was input with no internal generation. If data is generated 
internally, the default printing is all input data cards and 

a second printing of input data plus generated data. 

(4) The boundary condition code is used to designate those nodes which 
will have fixed values of temperature in the solution process. The 
fixed value of temperature is entered in the T(N) array. 

(5) IVhen CT (Col. 1) is equal to the character "C," the values input 
in columns 36-65 are interpreted as the cylindrical (R,Y,0) 
coordinates of node "N." Y is the axis of symmetry. R is the 
distance of a point from the Y-axis. The angle 0 is measured 
clockwise from the positive Z-axis when looking in the positive 
Y direction. The cylindrical coordinate values are printed as 
entered on the card, but immediately after printing the global 
Cartesian values are computed from the input entries . 

(6) Nodal point cards need not be input in node-order sequence; 
eventually, however, all nodes in the integer set (1, NUMNP} 
must be defined. Nodal data for a series of nodes 
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III. NODAL POINT DATA (concluded) 


{Ni, Ni + (1 X KN2), Ni + (2 X KN2), . . N2> 


may be generated from information given on two cards in sequence: 


CARD 1 / Ni, ID(Ni), X(Ni), . . KNi, T(Ni) / 

CARD 2 / N2, (D(N2), XCN2), . . KN2, TCN23 / 


KN2 is the mesh generation parameter given on the second card of 
a sequence. The first generated node is Nj + (1 x KN2) ; the 
second generated node is Ni + (2 x KN2) , etc. Generation 
continues until node nmber N2 - KN2 is established. Note that 
the node difference N2 - Ni must be evenly divisible by KN2 • 
Intermediate nodes between Ni and N2 are' located at equal 
intervals along the straight line between the two points. Boundary 
condition codes for the generated data are set equal to the values 
given on the first card. Node temperatures are found by linear 
interpolation between T(Ni) and T(N£) . Coordinate generation 
is always performed in the (X,Y,Z) system, and no generation is 
performed if KN2 is zero (blank) . 


NI N2 

e • e e e « 9 a 0 


(7) The nodal temperatures are used to specify temperatures which are 
fixed in the solution process. In the nonlinear analysis, the 
first iteration is performed with the thermal parameters evaluated 
for all nodes at a zero temperature value. Thereafter nodal 
temperatures including the fixed values specified here are used 
to evaluate the thermal parameters for the next iteration. 
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IV. ELEMENT DATA 


TYPE 1 - CONDUCTION/CQNVECTION ROD ELEMENT 

Rod elements Cfig- 5) are identified by the number 1. A linear 
temperature variation is assumed between nodes. Internal heat 
generation, prescribed surface heating or convective surface heating 
are incorporated into the element formulation. 

A. Control Card (815, 8A5) 

Columns 1-5 The number 1 

6-10 Total number of rod elements in this element group 
11-15 Number of material property cards 
16 - 40 Blank 

41 - 80 Any desired identification to be printed with output 

B. Material Property Cards (2I5,F10.0) 

The program expects the number of material property cards given 
above. 

Notes Columns 1-5 Material identification number 
(1) 6-10 Table number for thermal conductivity 

(nonlinear analysis only) 

11 - 20 Thermal conductivity (k) (required for 
linear analysis only) 

Surface heat fluxes (positive into surface) (BTU/hr) . 


X (local axis) 


Positive sign convention for conduction 
heat flux (BTU/hr) 


C. Element Data Cards 

One card per element in sequential order of element number starting 
with one. If there is surface heating or convection heat transfer 
two cards are required. 
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IV. ELEMENT DATA (continued) 


Card 1 (Required) (SI5,3F10.0) 


Notes 

Columns 1 - 

5 

Element number 

(2) 

6 - 

10 

Node number I 


11 - 

15 

Node number J 


16 - 

20 

Material identification number 

(3) 

21 - 

25 

Optional element generation parameter KG 
for automatic generation of element data 


26 - 

35 

Cross-sectional area for conduction 


36 - 

45 

Heat generation per unit volume (e.g., 
BTU/HR-FT-^*3) 

(4) 

46 - 

55 

Area factor for surface heating or convection 


Card 2 (Optional - required onl/ if the area factor is greater than 
zero) (5F10.0) 


Notes 

Columns 1 

- 10 


11 

- 20 


21 

- 30 

(5) 

31 

- 40 


41 

- 50 


Specified surface heat transfer rate 
(e.g., BTU/HR-FT*’'2) (positive into element) 
Convective medium heat transfer coefficient 
Hj at node I 

Convective medium temperature Tj at node I 
Convective medium heat transfer coefficient 
Hj at node J 

Convective medium temperature Tj at node J 


NOTES/ 

(1) For a linear analysis the thermal conductivity k input 

on the material property card is used to compute the thermal 
conductance matrix and the heat flux recovery matrix for an 
element. For a nonlinear analysis and a table number greater 
than zero, the thermal conductance matrices are initially com- 
puted using k as unity. Later, after the thermal parameter 
tables have been read in, the matrices are multiplied by 
appropriate values of k determined from the parameter tables. 

The temperature used in the table is the average temperature 
of the element, i.e., (Tj + Tj)/2 . 

(2) The order of I and J determines the direction of the local 
X-axis (see fig. 5), Conduction heat fluxes are positive in 
the direction of the local X-axis. 

(3) If a series of elements exists such that the element number, , 
is one greater than the previous element number (i.e., 

Nj^ = Nj^_i T 1) and the nodal point number can be given by 

li = Ii_i - KG 
Ji = Ji_i - KG , 
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IV. ELEMENT DATA (continued) 


then only the first and last elements in the series need be provided. 
The material identification number and the temperature for the 
generated elements are set equal to the values on the last card. 

If KG (given on the last card) is input as zero, it is set to 
one by the program. 

(4) If the area factor is greater than zero, the second card will be 
read. The area factor is used to compute the surface area for 
surface heat transfer, i.e., A (surface) = Area factor length 
of element. 

(5) If Hj or Tj is left blank, the program will set 
Hj = Hi and Tj = Ti . 



SURFACE HEATING 



Figure 5. Conduction/convection rod element. 
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IV. 

TYPE 


ELEMENT DATA (continued) 


3 - CONDUCTION/CONVECTION QUADRILATERAL ELEMENT 

Quadrilateral elements (fig. 6) are identified by the ntunber 3. The 
element is based on an isoparametric formulation. The nodes can be 
located at general points in space, but they must lie in a plane. The 
element conduction heat fluxes are computed at the element centroid 
in local coordinates . The element may be laminated with an arbitrary 
number of different layers with different conduction tensors for each 
layer. Internal heat generation, prescribed edge or surface heating, 
or convective heating on all four edges and the top and bottom surfaces 
of the element are included in the element. 


A. Control Card (815, 8A5) 


Columns 1-5 
6-10 
11 - 15 
16 - 40 
41 - 80 


The number 3 

Total number of quadrilateral elements in this group 

Number of material property card sets 

Blank 

Any desired identification to be printed with output 


B. Material Property Card Sets 


One card is required per set. Cards 2, 3, 
optional . 


number of laminae is 


Card 1 

(3I5,5F10 

.0) 


Notes 

Columns 

1 - 

5 

CD 


6 - 

10 

(2) 


11 - 

15 



16 - 

25 

(3) 


26 - 

35 



36 - 

45 



46 - 

55 



56 - 

65 


Material identification number 
Table number for thermal conductivity 
tensor temperature variation (nonlinear 
analysis only) 

Number of laminae 
Lamina thickness, t 
Conductivity tensor component, 

Kxx 

Conductivity tensor component. 


^xy 

Conductivity tensor component. 


< Lamina 1 


K- 


■YY 


Material axis angle, 0 (degrees) / 
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IV. ELEMENT DATA (continued) 


Conduction fluxes 



Card (2j 3, . . . No. of Laminae) (ISXj SFiO.O) 


Columns 16 - 25 
26 - 35 

36 - 45 

46 - 55 

S6 - 65 


Lamina thickness 
Conductivity tensor component , 

Conductivity tensor componentj 
^xy 

Conductivity tensor component^ 
^^^erial axis angle, 0 


Lamina 2, 3, . . , 
No . o£ Lafliinae 


C. Element Data Card Sets 


One card per element is required in increasing numerical order. 
Missing elements are generated. If there is edge or surface heating, 
additional element cards are required. 
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IV. ELEMENT DATA Ccontinued) 


NOTES/ 


Card 1 (9I5,F10.0) 


Notes 

Coltunns 1 

- 5 

Element number 

(4) 

6 

- 10 

Node I 


11 

- 15 

Node J 


16 

- 20 

Node K 


21 

- 25 

Node L 


26 

- 30 

Material identification number 


31 

- 35 

Element generation parameter KG 


36 

- 40 

lEDGE .EQ.O No edge heating 


ox edge convection 

.EQ.1,2,3,4 Number of edges for 
which there is edge 
heating or edge 
convection 

41 - 45 ISURF .EQ.O No surface heating or surface 

convection 

.EQ.l Heating or convection on top 
surface 

.EQ.2 Heating or convection on top 
and bottom surfaces 

46 - 55 Volumetric heat generation rate (e.g., 
BTU/HR-FT^) 

Card Set 2 [lEDGE cards) (2I5,5F10.0) 


Notes 


( 6 ) 


Columns 1-5 
6-10 
11 - 20 

21 - 30 
31 - 40 
41 - SO 
51 - 60 


Edge node, N1 
Edge node, N2 

Edge heat loading, q, (©-g-j BTU/HR-FT^) 

(heat flux is positive into element) 
Convection coefficient, hi, at node N1 
Convective medium temperature, Tj^, at node N1 
Convection coefficient, h. 2 , at node N2 
Convective medium temperature, Tz, at node N2 


Card Set 3 (ISURF cards) (8F10.0) 


Notes 

Columns 1 

- 10 

Convection 

(7) 



convective 


11 

- 20 

Convective 


21 

- 30 

Convection 


31 

- 40 

Convective 


41 

- 50 

Convection 


51 

- 60 

Convective 


61 

- 70 

Convection 


71 

- 80 

Convective 


coefficient hj at node I, or 

surface heating, 

medium temperature Tj at node I 
coefficient, hj at node J 
medium temperature Tj at node J 
coefficient, hj^, at node K 
medium temperature Tj^ at node K 
coefficient, hL, at node L 
medium temperature, Tj;^, ar node L 


(1) All of the components of the conductivity tensor are assumed to have 
the same temperature variation in a nonlinear analysis so that only 
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IV. ELEMENT DATA Ccontinued) 


one table is input for the entire tensor. The look-up temperature 
is (Tj + Tj + Tfj + Tl)/ 4 . For a nonlinear analysis, the element 
conductance matrix is formed for the first iteration using the 
conductivity tensor entered as input data. On subsequent 
iterations the thermal conductivity table is used as a multiplier 
of this tensor. Thus, for a single layer the user may 
input a conductivity tensor with the largest value normalized to 
1.0 and enter the actual conductivity values in the table. 

(2) For an element with one homogeneous layer, only the first card is 
required . 

(3) For an isotropic material, the conductivity value K should be 

entered as . The remaining entries may be left blank. The 

program will set K^y. = 0 , Kyy = . 

(4) The orientation of the local X-axis is from I to J (see fig. 6) . 

The local y-axis then lies in the IJKL plane, and the direction of 
the local z-axis is determined by the right-hand rule. Element 
conduction heat fluxes are positive in the local coordinate system. 

(53 Element cards must be in element number sequence. If cards are 
omitted, element data will be generated. The node numbers will 
be generated with respect to the first card in the series as 
follows. 


^n = Vl 


KG 

Jn = J^n-l 


KG 

^n “ ^-1 


KG 

^ ~ ^n-1 


KG 


All other element information will be set equal to information on 
the last card. 

» 

(6) If h .2 and T 2 are left blank, the program will set h 2 = hj , 

T2 = Ti . 
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IV. ELEiVENT DATA (continued) 

(7) The top surface is located on the positive local z axis: 


r 



Surface heating is positive into the element. If there is surface 
heating, the heat flow is entered in columns 1-10, and the remainder 
of the card is blank. For uniform convection, if hj , Tj , etc. 
are left blank, the program will set 

hj = hK = hL = hi 
Tj = Tk = Tl = Tl . 
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Figure 6. Conduction/ convection quadrilateral element. 
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IV. ELEMENT DATA (continued) 
TYPE 8 - MASS -TRANS PORT ELEMENT 


Mass transport elements (£ig. 7) are identified by the number 8. The 
element is used to represent convective energy transport due to a mass 
flow rate m . 


A. Control Card (815, 8A5) 


Columns 1 - 5 

6 -'10 
11 - 15 
41 - 80 


The number 8 

Total number of elements in this group 
Number of thermal -fluid property card sets 
Any desired identification to be printed with the 
element output data 


B. Fluid Properties (I5,F10. 0, 15) 


Notes 

CD 


Columns 1 - 5 

6-15 
16 - 20 


Property identification number 

Fluid specific heat, Cp 

Table number for fluid^specific heat 


C. Element Data Cards (5I5,F10.0) 


One card per element is required in increasing numerical order. 
Missing elements are generated. 


Notes 

Columns 1 - 

5 

Element number 

(2) 

6 - 

10 

Node number, I 


11 - 

15 

Node number, J 


16 - 

20 

Property identification number 

(3) 

21 - 

25 

Element generation parameter, KG 


26 - 

35 

Fluid mass flow rate (e.g., Ibm/hr) 



NOTES/ 

(1) For a linear analysis, the fluid specific heat read-in on the 

fluid property card is used in element computations. For a non- 
linear analysis, values from a specific heat table are used if 
the table number is greater than zero. 
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IV. ELEMENT DATA (continued) 


(2) The nodal coordinates are arbitrary and are used only in the plot 
output. The order o£ the element nodes determines the direction 
of fluid flow, i.e., the fluid flow is from node I to J. 

(3) Missing elements are generated using the same scheme as for the 
rod element, i.e., node numbers will be generated with respect 
to the first card as follows: 


i^ ='ii-i + kg 
>^1 = Ji-l + kg . 


All other element information will be set equal to data from the 
last card. 
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LINEAR FLUID 



Figure 7. Mass transport element. 
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IV. ELEMENT DATA (continued) 

TYPE 9 - SURFACE-CONVECTION ELEMENTS 


Surface-convectxon elements (fig. 8) are identified by the number 9. The 
elements are used to represent convective heat transfer between a surface 
and a fluid with unknmvn temperature. Two elements, a quadrilateral 
and a triangle, are available. 


A. Control Card (815, 8A5) 


Columns 1 - S 
6-10 
I 

11 - 15 
41 - 80 


The number 9 

The number of surface convection elements in this 
group 

Number of thermal -fluid property card sets 
Any desired identification to be printed with 
the element output data 


B. Fluid Properties (I5,F10 .0,15) 


Notes 

( 1 ) 


Columns 1 - 5 

6-15 
16 - 20 


Property identification number 

Convection coefficient, h 

Table number for convection coefficient 


C. Element Parameter Data Cards (7I5,F10.0) 

One card per element is required in increasing numerical order. 
Missing elements are generated. 


NOTES/ 


(13 


Notes 

Columns 1 - 

5 

Element number 

(2) 

6 - 

10 

Node number, I 


11 - 

15 

Node number, J 


16 - 

20 

Node number, K 

(3) 

21 - 

25 

Node number, L (default .EQ.O) 


26 - 

30 

Property identification number 

(4) 

31 - 

35 

Element generation parameter, KG 

(5) 

36 - 

45 

Area factor for convection (default .EQ.l 

For 

a linear analysis. 

the fluid convection coefficient read- in on 

the 

fluid property 

card 

is used in element computations. For a 


nonlinear analysis, values from a convection coefficient table are 
used if the table number is greater than zero. 


(2) For the quadrilateral element, nodes I and J always denote fluid 
nodes; for the triangle, I denotes the fluid node, 

(3) For a triangular element, leave L blank or enter L as zero. 

(4) Missing elements are generated using this same scheme as for the 
quadrilateral conduction element, i.e., node numbers will be 
generated with respect to the first card as follows: 
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IV. ELEMENT DATA (continued) 


II = Ii_i + KG 

Ji = Ji_i + KG 

Ki = + KG 

Lj_ = -i- KG . 

All other element information will be set equal to data from the 
last card. 

(5) The area factor is used to compute the convection surface area, 

e.g., for the quadrilateral, A (surface) = area factor * distance K-L 
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CONVECTION SURFACE 



TYPICAL SURFACE NODE 


TYPICAL FLUID NODE 


(a) Quadrilateral Element. 



(b) Triangular Element. 


Figure 8. Surface convection elements with unknown 
fluid temperatures. 
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IV. ELEMENT DATA (continued) 

TYPE 10 - TUBE/FLUID INTEGRATED ELEMENT 


Tube/fluid integrated elements (fig. 9) are identified by the number 10. 
The element represents conduction/convection heat transfer in a thin tube 
of constant thickness and flow area enclosing a fluid with mass flow 
rate m . Heat loading on the tube external surface due to a specified 
heating or a convective exchange with a surrounding media is included. 
Pressure drop computations are performed as an option. 

A. Control Card (815, 8A5) 


Columns 1 
6 
11 
16 


21 

41 


5 The number 10 

10 Total number of tube/fluid elements in this group 
15 Number of thermal-fluid property card sets 
20 Flag for pressure drop calculations 

.EQ.O; Pressures are not calculated 

.GT.O; Pressures are calculated 

40 Blank 

80 Any desired identification to be printed with the 
element output data 


B, Thermal-Fluid Property Card Sets 


Card 1 

- Tube Properties 

(I5,F10.0,I5) 

Notes 

Columns 1 

- 5 

Property identification number 

CD 

6 

- 15 

Thermal conductivity, k 


16 

- 20 

Table number for tube thermal conductivity 

Card 2 

- Fluid Properties 

(F10.0,2(F10.0,I5),I5) 

Notes 

(2) 

Columns 1 

- 10 

Fluid convection coefficient, h 

C3) 

11 

- 20 

Ejtponent in equation for modification of 


21 

- 25 

convection coefficient, n (real) 

Table number for convection coefficient 

C4) 

26 

- 35 

Fluid specific heat, Cp 


36 

- 40 

Table number for fluid specific heat 

C5) 

41 

- 45 

Table number for fluid viscosity 


(Optional) Card 3 - Tube-Fluid Properties for Pressure Recovery 





(2F10, 

,0,2(F10.0,I5),2F10.0) 

Notes 





(6) 

Columns 

1 

- 10 

Tube hydraulic diameter, D^ 

(7) 


11 

- 20 

Fluid friction factor, f 

(8) 


21 

- 30 

Exponent in equation for correction of 
friction factor, m (real) 



31 

- 35 

Table number for fluid friction factor 
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IV. ELEMENT DATA (continued) 


Notes 


( 9 ) 

Columns 36 

- 45 

Fluid density, p 


46 

- 50 

Table number for fluid density 


51 

- 60 

Gas constant, R 


61 

- 70 

Proportionaility constant in Newton's 
second law, g^ 


C. Element Data Cards 


One card per element is required in increasing numerical order. 
Missing elements are generated. If there is external surface heating 
on the tube, two cards per element are required. 


Card 1 - Element Parameters (8I5,4F10.0) 


Notes 

Columns 1 

- 5 

Element number 


(10) 

6 

11 

- 10 
- 15 

Node number, I i 
Node number, J / 

Fluid nodes 


16 

21 

- 20 
- 25 

Node number, K l 
Node number, L > 

Tube nodes 


26 

- 30 

Property identification number 

(11) 

31 

- 35 

Element generation parameter, KG 

(12) 

36 

- 40 

ISURF.EQ.O; 

No surface heating or 
convection 




,GT.0; 

Surface heating or 
convection 


41 

- 50 

Tube cross-sectional conduction area 

(13) 

51 

- 60 

Perimeter of tube for internal convective 




heat transfer to 

fluid 


61 

- 70 

Fluid mass flow 

rate (e.g. , Ibm/hr) 

(14) 

71 

- 80 

Element inlet pressure, Pj 


Card 2 - Exteimal Tube Heating or Convection Data (6F10.0) 
Notes 


(15) 

Columns 1 - 

10 

Area factor for surface heating or 
convection (default .EQ.1.0) 


11 - 

20 

Specified surface heating rate (e.g., 
BTU/HR-FT‘'''*2) (positive into surface) 


21 - 

30 

Convective heat transfer coefficient, hj^ 
at node L 


31 - 

40 

Surrounding medium temperature, Tl^ at 
node L 

(16) 

41 - 

50 

Convective heat transfer coefficient, hj^ 
at node K 


51 - 

60 

Surrounding medium temperature, Tj^, at 
node K 
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IV. ELEMENT DATA (continued) 
NOTES/ 


(1) The thermal conductivity is used to represent the axial conduction 
of heat in the tube wall. The thermal conductivity of the tube 
wall may be constant or may be entered in tabular form for a 
nonlinear analysis. The look-up temperature is (Tj^ + Tl)/2 . 

(2) The fluid convection coefficient h is used to represent 
convective heat transfer between the tube and fluid. The 
convection coefficient may be constant or may be entered in 
tabular form for a nonlinear analysis . The look-up temperature 
IS (Tj + Tj)/2 . 


(3) In the nonlinear analysis, the convection coefficient may be 

modified at each iteration for a variation of fluid temperature 
at the flow section. The correction takes the form: 


Gases : 


h> = h (Tb) 



n 


Liquids : 


h' = h (Tb) 



n 


For a gas, T^^ denotes the wall temperature and Tb denotes the 
bulk fluid temperature. For a liquid, denotes the viscosity 

evaluated at the wall temperature and Ub denotes the viscosity 
evaluated at the fluid bulk temperature. No modification is 
performed if the exponent n (real) is entered as blank or zero. 

(4) The specific heat C_ is used in representing the heat transfer 
due to mass transport and may be entered as a constant or as a 
tabular function of temperature. 

(5) The fluid viscosity is required only if the correction described 
above in note (3) is to be performed for a liquid. Otherwise, a 
table number for viscosity is not required. 

(6) The tube hydraulic diameter is defined by 


Dji = 4 Flow cross-sectional area 
wetted perimeter 


(7) The fluid friction factor f is used to compute the pressure drop 
in an element. The pressure drop is computed from the equation. 
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IV. ELEMENT DATA (continued) 


AP = £ 


g2 


1 ^ G2 


Pm Sr 



where 


AP - pressure drop, (Pj - Pj) 
f - friction factor 
L - element length 

- hydraulic diameter 

G - mass flow rate/flow area (e.g., lbm/hr/£t^) 

Pjjj - element mean density, (pj + Pj!)/2 

pj, pj - fluid densities evaluated at the temperatures of 
the fluid at nodes I, J 

- proportionality constant in Newton's second law 

/ 32,17 ft - Ibm \ 

(^e.g., gc ) 


(8) In the calculation of pressures for a nonlinear analysis, the 

friction factor may be modified for a variation of fluid temperature 
at the flow cross section. The correction takes the form; 


Gases : 

11 

Liquids : 

f = f (Tb) 



in 


m 


where the quantities have the same definition as in note (3) . 

The modification is not performed if the exponent m (real) is 
entered as blank or zero. 

(9) Pressure drops are computed for three density cases; (1) constant 
density, (2) variable density as specified by a density-temperature 
table, and (3) an idea gas. If the density table number is entered 
as zero, case (1) is assumed. If the density table number is 
greater than zero, case (2) is assumed. 
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IV. ELEMENT DATA (continued) 


If the gas constant R is entered as a positive quantity, case (3) 
is assumed. For case (3) the pressure drop equation above is 
solved simultaneously with the gas law P = pRT . 

(10) The direction of fluid flow is determined by the node numbering 
sequence. Flow is from node I to node J (see fig. 9) . 

(11) Element cards must be in element number sequence. If cards are 
omitted, element data will be generated. The node numbers will 
be generated with respect to the first card in the series as 
follows ; 


^n 

= In-1 

-f- 

KG 

Jn 

" Jn-1 

+ 

KG 

Kn 

= Kn-i 

+ 

KG 

In 

® In-1 

+ 

KG . 


All other information will be set equal to the data on the last 
card, 

(12) ISURF.GT.O indicates the tube is heated externally by a specified 
heat flux or convectively. The program expects to read a second 
card with the heating data. 

(13) The perimeter of the tube is used to compute the wetted area for 
convective heat transfer to the internal fluid by multiplying it 
by the element length. 

(14) Pressures are computed at successive nodes by Pj = Pi - AP 
until a new inlet pressure is specified for an element. 

(15) The surface area for external heating is computed as the product 
of the area factor times the perimeter times the element length. 

(16) If h}j and Tj(; are left blank, the program will set ht^ = h^- 
and Tk = Tl . 


EEPRODUCIBILITY of the 

original page is poor 
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Figure 9. Integrated tube/£luid element. 
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IV. ELEf^IENT DATA Ccontinued) 

TYPE 11 - PLATE-FIN/FLUID INTEGRATED ELEME^^^ 


Plate-£in/ fluid integrated elements 10) are identified by the 

number 11. The elements represent conduction/convection heat transfer 
in a coolant passage defined by two plates connected by an internal 
fin. Fluid flows through the passage with mass flow rate m . 


Control 

Card (815 

,8A5) 

Columns 

1-5 

The number 11 


6-10 

Total number of plate-fin/fluid elements in this 
group 


11 - 15 

Number of thermal-fluid property card sets 


16 - 20 

Flag for pressure calculations 

.EQ.O] Pressures are not calculated 

.GT.O; Pressures are calculated 


21 - 40 

Blank 


41 - 80 

Any desired identification to be printed with the 
element input data echo. 


B. Thermal -fluid Property Card Sets 


Card 1 

- Fin Properties (I5,F10.0,I5) 

Notes 

Columns 

1-5 

Property identification number (Fin and 
fluid properties) 

CD 


6-15 

Thermal conductivity 



16 - 20 

Table number for fin thermal conductivity 

Card 2 

- Fluid Properties 

(F10.0,2(F10.0,I5) ,15) 

Notes 




C2) 

Columns 

1-10 

Fluid convection coefficient, h 

(3) 


11 - 20 

Exponent in equation for modification of 
convection coefficient, n (real) 



21 - 25 

Table number for convection coefficient 

C4) 


26 - 35 

Fluid specified heat, Cp 



36 - 40 

Table number for fluid specific heat 

(5) 


41 - 45 

Table number for fluid viscosity 

(Optional) Card 

3 - Properties for Pressure Calculations 



(2F10.0,2(F10.0,I5),2F10.0) 

Notes 




C6) 

Columns 

1-10 

Hydraulic diameter, Dj| 

(7) 


11 - 20 

Fluid friction factor, f 

(8) 

• 

21 - 30 

Exponent in equation for modification of 
friction factor, m (real) 



31 - 35 

Table number for fluid friction factor 
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IV. ELEMENT DATA (continued) 


Notes 

(9) 


Columns 36 - 45 
46 - 50 
51 - 60 
61 - 70 


Fluid density, p 

Table 'number for fluid density 

Gas constant R 

Proportionality constant in Newton's 
second law, 


C. Element Parameter Data Cards 


Two cards per element are required in increasing numerical order. 
Missing elements are generated. 

Element Parameters (10I5,2F10.0,/,3F10.0) 

Card 1 


Notes 

Columns 1 - 

5 

(10) 

6 - 

10 


11 - 

15 


16 - 

20 


21 - 

25 


26 - 

30 


31 - 

35 


36 - 

40 

(11) 

41 - 

45 

(12) 

46 - 

50 



51 

- 60 

(13) 

61 

- 70 

Card 2 



Notes 

(14) 

Columns 1 

- 10 

(15) 

11 

- 20 


21 

- 30 

(16) 

31 

- 40 


Element nimber 
Node number , I 
Node number, J 
Node number, K (Fluid node) 

Node number, L 
Node number, M 

Node number, N (Fluid node, inlet) 

Property identification mmiber 

Element generation parameter, KG 

Flag for fin efficiency 

.EQ.O; Fin efficiency computed 

.NE.Oj Fin efficiency set equal to one 

Fluid mass flow rate (e.g., Ibm/hr) 

Element inlet pressure, Pj^j 


Effective fin thickness 

Effective width of top wall for convection 
(default, 1.0) 

Effective width of bottom wall for 
convection (default, 1.0) 

Fin area factor (default, 1.0) 


NOTES/ 

(1) The thermal conductivity is used to calculate two-dimensional 

heat conduction in the fin. The fin connects the top and bottom 
walls, and heat conduction is represented by an isoparametric 
quadrilateral finite element formulation. The thermal conductivity 
may be constant or entered in tabular form for a nonlinear analysis . 
The look-up temperature is (Tj + Tj + T^ + T]V()/4 . 
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IV. ELEMENT DATA Ccontinued) 


(2) The fluid convection coefficient h is used to represent convective 
heat transfer between the top and bottom walls and between both 
sides of the fin and the fluid. The convection coefficient may be 
constant or be entered in tabular form for a nonlinear analysis . 

The look-up temperature is (Tj^j + Tj() /2 . 


(3) In the nonlinear analysis ^ the convection coefficient may be 

modified at each iteration for a variation of fluid temperature at 
the flow section. The correction takes the form: 


Gases : 


h' = h (Tb) 



n 


Liquids: h' = h (T^) 



For a gas, Tw denotes the wall temperature and Tjj denotes the 
bulk fluid temperature; for a liquid, denotes the viscosity 

evaluated at the wall temperature and y^ denotes the viscosity 
evaluated at the fluid bulk ten^jerature. No modification is 
performed if the exponent n (real) is entered as blank or zero. 

(4) The specific heat Cp is used in representing the heat transfer 
due to fluid floii? anu may be entered as a constant or as a tabular 
function of temperature. 

(5) The fluid viscosity is required only if the modification described 
above in note (3) is to be performed for a liquid. Otherwise, a table 
number for viscosity is not required. 

(6) The passage hydraulic diameter is defined by 


Dh = 4 


Flow cross-sectional area 
wetted perimeter 


(7) The fluid friction factor f is used in computing the pressure 
drop in an element. The pressure drop is computed from the 
equation 


AP 


Jl ^ _L + ^ 

% 2gc Pm Sc 




EEPEODUCIBILITy OF THE 
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IV. ELEMENT DATA (continued) 
where 


AP 

£ 

L 

% 

G 

Pm 
PK> PN 

Sc 


pressure drop, (Pj^ - pj() 
friction factor 
element length 
hydraulic diameter 

mass flow rate/flow area (e.g. , Ibm/hr/ft^) 
element mean density, (pj/ + Pjyj)/2 

fluid densities evaluated at the temperatures of the 
fluid nodes K, N 

proportionality constant in Newton's second law 


( 


e.g. , 


®c 


32.17 ft - Ibm 

Ibf - sec2 


(8) In the calculation of pressures for a nonlinear analysis, the 
friction factor may be corrected for a variation of fluid 
temperature at the flow cross section. The correction takes the 
form: 


Gases : 


Liquids : 


where the quantities have the same definition as in note (3). The 
modification is not performed if the exponent m (real) is 
entered as blank or zero. 

(9) Pressure drops are computed for three density cases. (1) constant 
density, (2) variable density as specified by a density-temperature 
table, and (3) an ideal gas. If the density table number is 
entered as zero, case (1) is assumed. If the density table number 
IS greater than zero, case (2) is assumed. 

If the gas constant R is entered as a positive quantity case (3) 
IS assumed. For case (3) the pressure drop equation above is 
solved simultaneously with the gas law P = pRT . 




f = £ (Tb) 



1 Pw 

f = f (Tb) 



m 


EEPRODUCIBILITY OF THE 
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IV. ELEMENT DATA (continued) 


(10) The direction of fluid flow is determined by the node numbering 
sequence. Flow is from node N to node K (see fig. 10). 

(11) Element cards must be in element number sequence. If cards are 
omitted, element data will be generated. The node numbers will be 
generated with respect to the first card in the series as follows: 


In = In-l + kg 
■In ~ ’In-1 kg 

Kn = Kn-i + KG 
Ln = Ln-1 + KG 
Mn = Mn-i + KG 
Nn = Nn-i + KG 


All other information will be set equal to the data on the last 
card. 

(12) The fin efficiency n is computed by the equation 


2 cosh m& - 1 
mJ?- sinh m£ 


where I is the average height of the fin and 



The fin efficiency is used to modify the convective heat transfer 
between the fin and fluid for the linear temperature distribution 
assumed in the surface convection finite element (see note (16) 
below) . 

(13) Pressures are computed at successive nodes by Pk = ?N " 
until a ne^7 inlet pressure is specified for an element. 

(14) The fin thickness is used in two ways. The thickness is used in 
representing the conduction heat transfer of the fin. In addition, 
the fin thickness is subtracted from the widths of the top and 
bottom walls in the computation of convection areas. For multiple 
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IV. ELEMENT DATA (concluded) 


fins, an effective fin thickness equal to the number of fins 
times the thickness of a single fin should be used. 

(15) The top and bottom widths are used to compute the convection areas 
from the walls to the fluid (see note 14) . The average of these 
widths is also used in the computation of the flow areas at the 
inlet and outlet of an element. 

(16) The fin area factor for convection may be used to account for 
multiple fins. The fin surface area is multiplied by this 
factor. The convective heat transfer between the fin and fluid is 
based upon the equation 


q = r|h (A3 * Ap) 


■ Ti Tj + Tl + Tm 

4 '^BULK 


where 

ri = fin efficiency (see note 12) 
h = convection coefficient 
Ag = surface area of fin (2 sides) 
Ap = fin area factor 
'I’bULK = (Tn + Tk)/2 . 
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(a) Side View. 





FLOW 


(b] Oblique View. 


Figure 10. Integrated plate-fin/f luid element. 
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V. THERMAL PARAMETER TABLES 

Thermal parameter tables are required for nonlinear thermal analysis. The 
total number of thermal parameter tables is entered on the master control 
card as NUMTB (see section II) . Individual table numbers for reference 
to the data input here are read in as part of the element input data. 

The thermal parameter data tables are described by the following sequence 
of data cards : 


A. Control Card (2I5,7A10) (one card for each table) 

Columns 1-5 Table number 

6-10 Number of data points given in table 
11 - 80 Any desired table heading information 

B. Thermal Parameter Table (SFIO.O) (4 points per card, as many cards as 
required) (typical card) 


Columns 1-10 
11 - 20 
21 - 30 
31 - 40 
41 - 50 
51 - 60 
61 - 70 
71 - 80 


Temperature for point 
Thermal parameter for 
Temperature for point 
Thermal parameter for 
Temperature for point 
Thermal parameter for 
Temperature for point 
Thermal parameter for 


1 

point 


Point 

2 

point 

2 ) 

Point 

3 

point 

3I 

Point 

4 

point 

4 ) 

Point 
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APPENDIX C 


INPUT DATA FOR TAPPLT 
General Setup of Deck 

In general the input deck for the plotting program consists of three 
separate groups of data as sho\m schematically in figure 11. These groups 
are as follows : 

(1) A single card containing any desired title information, 

(2) NA^iELIST OPTION containing values to allocate storage in blank 
common and control values specifying various program options, and 

(3) NAiMELIST PICT containing values to specify the type of plot 
desired and information to be included on the plots . 

Input Data Cards 

I. HEADING CARD - This single card contains any desired alphanumeric 

information in columns 1 to 80. The title will appear at the beginnin 

I 

of the plots. 

II. NAMELIST OPTION - This NAMELIST contains values to allocate storage in 
blank common and control values specifying various program options. 


FORTRAN 

name 

Default 

value 

Description 

NNDEST 

200 

Estimated number of nodes, must be equal 
to or greater than the actual number of nodes 

NWDISP 

0 

0 no temperature data 

1 temperature data 

KPLOT 

1 

Specifies the type of output device to be used 

1 CalComp 

2 CalComp with plotting speed reduced for 
Leroy pens 

3 VARIAN 

XSPACE 

10.0 

Space between plots in x~direction, in inches 


E^iPKODUCIBILm OF 'lUi^ 
PAGE K POOR 




III. 


FORTRAN Default 

name value Description 


PSIZE 25.0 Paper size in y-direction, in inches [used in 

scaling of plots to insure this dimension is 
not exceeded) 


NA^ffiLIST PICT - This NAMELIST contains values to specify the type of 
plots desired and the information that is to be included on the 
plots. 

FORTRAN 

name 

Default 

value 

Description 

KHORZ 

1 

Integer designating the horizontal axis of the 
viewing plane where 1 = Xq j ^ , and 

3 = Zo 

KVERT 

2 

Integer designating the vertical axis of the 
viewing plane where 1 = Xq » 2 = Yq , and 
3 = Zo 

PHI 

0.0 

Angular rotation of model about its x-axis in 
degrees [must be performed third) 

THETA 

0.0 

Angular rotation of model about its y-axis in 
degrees [must be performed second) 

PS I 

0.0 

Angular rotation of model about its z-axis in 
degrees [must be performed first) 

NEWFR 

1 

1 frame change before plotting [a frame change 

resets the x-origin past previous plot by 
XSPACE given in NAMELIST OPTION and resets 
the y-origin at 0.0) 

0 no frame change before plotting 

ISCALE 

1 

1 automatic computation of proper origin 
location and scaling of plot 

2 user-specified origin and scaling 

PLOTSZ 

10.0 

Maximum dimension desired on completed plot, in 
inches [used for scaling if ISCALE =1) 

XORGN 

0.0 

x-location of plot origin [used if ISCALE = 2) 

YORGN 

0.0 

y- location of plot origin [used if ISCALE =2) 

PSCALE 

■ 1.0 

Model size reduction factor [i.e., PSCALE is equal 


to actual model size divided by desired plot size 
(used if ISCALE = 2)) 


S^MODUCIBILITY OF THD 
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FORTRAN 

Default 



name 

value 


Description 

NOTAT 

0 

0 

no numbering on plots 



1 

numbering of nodes 



2 

numbering of elements 


XLHT 

0.15 

Height of integers specified by NOTAT, in inches 
(must be _> 0.07) 

KDISP 

0 

0 plot of thermal model 

1 plot of temperature surface 

2 exploded plot of model 

3 temperatures represented by vectors 

IDMAG 

2 

1 direct magnification of temperature data by 
DMAGS 

2 scaling of temperature data to a maximum 
value of DMAGS 

DMAGS 

1.0 

Magnification of temperatures (if KDISP = 1 or 3) 
Reduction factor of elements (if KDISP = 2) 

KSYMXY 

0 

1 symmetry about X-Y plane 

KSYMXZ 

0 

1 symmetry about X-Z plane 

KSYMYZ 

0 

1 symmetry about Y-Z plane 


Symmetries are performed consecutively (i.e., a plate quadrant with 
KSYMXZ and KSYMYZ equal to one would yield a complete plate) . 


XXMAXj 1.0 E+20 Locate cutting planes parallel to principal planes 

YYMAX, CX-Y, X-Z, Y-Z) to limit plot 

ZZMAX 


XXMIN, -1.0 E+20 
YYMIN, 

ZZMIN 


NDMAX 9999999999 Maximum node identification number to be 

included in plot 


NDMIN 0 Minimum node identification number to be 

included in plot 

NELMAX 9999999999 Maximum element identification number to be 

included in plot 


NELMIN 


0 Minimum element identification number to be 

included in plot 
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FORTRAN 

name 

Default 

value 


Description 

KODE 

0 

Specifies control option after plot is complete 

0 last plot;, exit from program 

1 read another NAMELIST PICT 


This section describes a complete basic set of input data if KODE = 0 
in NAMELIST PICT. 

The deck must end with NAMELIST PICT having a value of KODE = 0 . 
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APPENDIX D 


INPUT DATA AND PROGRAM OUTPUT FOR SAMPLE PROBLEMS 

Four sample problems are presented: fl) a linear conduction analysis 

of a transverse cross-section of a panel with a "D" tube, (2) a nonlinear 
analysis of a convectively heated, water cooled pipe, (3) a nonlinear analysis 
of a simplified heat exchanger, and (4] a linear analysis and plots of 
conduction in a simple fin. The sample problems are presented in figures 
12 to 15. Plotter output for sample problem (4) is sho™ in figure 16. 
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SAMPLE PROBLEM 1 


Linear Conduction Analysis of a D Tube Cross Section 
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INPUT DATA (SAMPLE PROBLEM 1) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


0 TUBE Ct^OSS SECTION 3 / 16/76 
58 2 0 1 0 0 


1 

0 

- 0.0825 

0.0 

0.5 

1 1 

0 

- 0.0825 

0.0 

0.25 

2 

0 

- 0.045 

0.0 

0.5 

12 

0 

- 0 . C 45 

0.0 

0.25 

13 

0 

- 0,0825 

0 . 0 

0.225 

1<4 

0 

- 0.045 

0.0 

0.225 

15 , 

0 

- 0.035 

0.0 

0.225 

16 

0 

0,0 

0.0 

0.225 

17 

0 

- 0.0825 

0.0 

0.190 

la 

0 

- 0.045 

0.0 

0.190 

19 

0 

- 0.035 

0.0 

0. 190 

20 

0 

0.0 

0.0 

0.190 

21 

0 

- 0.0825 

0.0 

0.160 

22 

0 

- 0.0450 

0.0 

0. 160 

23 

0 

- 0.035 

0.0 

0.160 

2 <V 

0 

0.0 

0.0 

0.160 

27 

0 

- 0.0825 

0.0 

0.120 

28 

0 

- 0.0450 

0.0 

0.120 

29 

0 

- 0.035 

0.0 

0.120 

30 

n 

0.0 

0.0 

0.120 

33 

o 

- 0.0825 , 

0 . 0 

0.080 
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PROGRAM OUTPUT (SAMPLE PROBLEM 1) 


.0. CROSS SECTION Jin 6>76._ 


Q. jN_L B_0 _L I . N. E ,0__R „K A T I ~a~ M 


NUMCER OF NOOAL POINTS 58 

NUnUER. OF .ELEMENT TYPES « ' 

NUMBER OF TABLES' o 

ANALYSIS LOOEINANA) ' * " 1 _ J L 

E0.0» DATA CHECK ONLY, 

EO.l, LINEAR, 

EQ.2, NONLINEAR ' ' ■'JJ" 

-PLOT COOE(NPLOT) « 'o 

EO.O, NO PLOTS generated" 

EO.l, UN0EFQPP6O PLOT 

. .EQ.2, TLMPERAIURE .PLOT 

ITERATION PAPAHETERS 

MAXIMUM ITERAHONS • 6 

- tolerance J,JLQfiA0EjL0CL 


BLANK COMMON LOCATIONS - 10967 


''a 

Is 3 


NODAL POINT INPUT DATA 


NOOr BOUNDARY CONDITION CODE 

NUMBER 

1 ----- ^ 

12 'T~ " " ■ 0 

13 0 

— 16 0 

15 0 

16 0 

- 17 - - 0 

la 0 

19 0 

- 20 _ 0 

21 0 

22 0 

- 23 0 

26 0 ■ ■ 

2 7 ~ 0 

- 2 a Q UlU 

20 . 0 

<0 0 

3 3 ^. 0 


nodal. B.0INT.j;00_RDINAIES „ 

^ Y . 7 kN TEMPERATURE 


r .033 

.. --.003 

=,06 5 

. _ r .065 

- -.033 

.= ,065 

-- -.035 

. 0.000 

-.08 3 

- -.065 . „ 

-.035 . 

_ 0.000 

-.063 . . 

. -.065 

- -. 035 . . 

.. 0.00 0 

.. -.03 3 

=-,^061 

- -.035 

0 . 000 
-.033 


0.000 

. 500 ,. — . 0 . 

0.0 00 

- 0.000 _ 

.250 „ 2 

o.ooc 

- 0,000 

-.500 0 

0.000 

. 0,000 

.250 2 . 

0,000 

- 0.000 . 

.225 0 

0.000 

- 0.000 

.225 0 

0.000 

0.000 

.225 -.0 . 

0.000 

. 0,000 

.225 0 .. 

0.000 

_ 0.000 

. 190 0 

0.000 

0.000 

.190 .. 0 

V * VyVJV / 

0.000 

0,000 

.190 0 

0.000 

.. 0,000 

..190 0 . _ 

0.000 - 

0.000 . 

.160 0 . 

0.000 

0.000 

.160 0 

0,000 

0.000 

.160 0 

0.000 

0.000 . . 

. 160 „ 0 

0.000 

. 0.000 

-. 120 . __ C 

0.000 . 

_ o . ooo . 

1 0 0 

0.000 

0.000 

.120 0 " 

C.JOO 

0 . 000 

* 1 ? L * 0 

0,000 

0.0 0 {, 

• 0 

o.coo 




34 

0 . 

-.04 5 

0.000. .. - 

,080_ __ 

0 . 

0,000 _ 


35 

0 

-.035 

0.000 

.080 

0 . 

0.000 


36 . 

0 

0.000 

0.000..__ 

.080 

0 . . 

0.000 


39 

0 

-.083 

0.000 

• 04 0._ 

0 

0.000 


40 

0 

-.045 

0.000 

,040._ 

0. 

0,000 . . . . 


41 

0 _ 

-.035 

0,000 

,040 

0 

0.00 0 ... 


47 

0 

-0.000 

0.000 

.040 

.0 

0.000 


45 

0 

-.083 

0.000 

0.000 

0 

O.OOO .. . 


46 

0 

r.045. 

0.000 

0.000 , 

0 

0.000 


47 

0 

-.035 

0.000 

O.OOO 

0 

0.000 . 


48 

0 . 

. —O.OOO. 

0.000 

.0.000 

o._ . 

0.000 

'c 

57 

0 

, .190 

0,000 

90.000 . . 

0. 

0,000 

c 

55 . 

0 „ . 

- .190 

0.000 

flO.OOO 

0 

0.000 . 

c 

53 

0 

. . . .190 

0.000 

. 70.000 .... 

0 „ ' „ 

0.000 

c 

51 

0 

. .190 

0,000 

_60.000 

0 . 

0.000 . . . 

c 

49 

0 . 

.19 0 

0.000 

_.50.000_ . 

0 

0.000 

c 

43 

0 „„ 

. ..190 ._ 

0.000 

„4 0.000 

0. _ 

0.000 

c 

37 

0 

.19Q._._ 

O.QQO- 

30. 000 

0 . 

0.000 

c_. 

31 

.o._ 

.190 

0_.000_.- 

70.000 

0 

0.000 


C_ 25 0 .190 0.000 10,000 0 0.000 

C _26 0 .225 0.000 10.000 0 0.000 

.C_.32 0 .22 5 Q.OOO 20.000 0 0.000 

C 36 0 .225 0.000. 30.000 0 0.000 

C .. 0 .22 5 0.000 40.000 0 0.000 „ 

,C 50 0 , 225 0,000 50.000 0 0,000 

C 52 0 .22 5 0.000 60.000 , . 0 0.000 

C 54 .. 0 . _ - . .22 5 0.000 70.000 0 0,000 

C 56 0 .225 0.000 80.000 0 0.000 

C 58 _ _ .0 .225 0.000 90.000 0 0.000 . 


generateo_nooal„oata 


.MODE ' _aOUND AR y'.CONO I T 1 0R_ ciToF Ni57L_tDIMjJxO~QRDT:MLT ES 


NUMBER 


X 


Z KN_ 

.tThperaiure 







1 

0 . 

-.083 

0.000 

.500 

0-000 

2 .. _.. 

0 ... 

-.0*5 

. 0.000. 

.500.__ 

0,000 „ 

3 . 


_-.083 

' 0.000 

_.,450 ._ _ 

. 0.000 . .. 

4 .. 

0 

-.045 

0.000 

.450, 

0.000 

5 

.0 

-.083 

. 0.000 . . 

. 4 OO 

0.000 

. .6 . . , 

. 0 

-.04 5. 

0.000 

.400 

0.000 

7 

0 

-.083 

0.000 

. 350 

0.000 

6 . 

0 ..... 

.. _-.045. 

0.000 

, 350 _ 

0.000 

9 

. . 0 . . 

. -.083 

0.000 

.300 

0.000 „ 

10 

0 

-.045 

0.000 

.300 

0.000 

_11 ._ . 

0 _.. 

-.083 

0.000 

_ .250 . 

0.000 

1? 

0 

“ < C 1 3 

0.000 

.250 

0.000 

13 

0 

-.0)3 

0.000 

.775 

0.000 


04 



1 A 

0 

-.045 

0.000 

.225 

_ 0 .._ 00 . 0 _. 

15 

0 

-.035 

0.000 

.22.5 

0.000 

16 

- 0 

0.000 

0.000 _ 

.2 25 ,._ 

0.000 

17 

Q 

-.083 

0.000 

.190 

Jl.OOO 

18 

0 

-.045 

0.000 

,190 

. 0.000 .. 

19 

0 

-.035 

0.000 

.190 

0 , 000 . . 

20 

0 

0.000 

0.000 

• J . 9.0 

. 0 . 000 ... 

21 

0 

-.083 

0.000 

_ . 160 

0.000 

22 

0 

-.045 

0.000 

.160 

o.oco 

23 

0 

-.035 

0.000 

.160 

0,00 0 .. 

2 ‘i 

0 

0.000 

0.000 

.160 _ 

0.000 

25 

0 

.033 

0.000 

.137 

0.000 _ 

2 5 

0 

.039 

0.000 

, 222 . 

0 . 000 . 

27 

0 

-.08 3 .,.. 

0.000 

. .120 

0,000 

26 

0 

-.045 

0.000 

.120 . _ . 

0.000 

29 

0 

-.035 

0.000 

.120 

0.000 

30 

0 

0 . 000 ... 

_ 0.000 

, 120 _ 

0,000 . 

31 

0 

.065 

. 0,000 

, 179 . 

0,000 

32 

0 

.077 

0 , 000 . 

,211 

0.000 

33 

0 

-.083 

0.000 

.080 . . .. 

0.000 

3 -i 

0 

-.045 

0.000 

.080 

0.000 

35 

0 

-.035 

0.000 

.080 

0,000 

35 

0 

0.000 

. 0.000 . _ 

.080 . . 

0.000 

37 

0 

.095 

0.000 

, 165 . . . . 

0,000 

38 

0 

'.112 

0,000 

195 

. 0 , 000 .. 

39 

0 

-.083 

0.000 

_ .040 . . 

0.000 

‘lO 

0 

-.045 

,. 0.000 . 

_., 040 . ... 

0.000 

A 1 

0 

-.035 

0.000 

.040 

0.000 

A 2 

0 

0.000 

. 0.000 

..040 

0,000 

43 

0 

. 122 

0.000 

, 146 . 

0.000 

44 

0 

.145 

0.000 

. 172 

0.000 . 

45 

0 

-.083 

0.000 

0.000 

0.000 

46 

0 

-.04 5 ..... 

0.000 

. 0.000 ... . 

0.000 

47 

0 

-.035 

o . ooo .... 

0.000 

.. 0.000 . 

48 

0 

o.noo 

. 0.000 

. . 0,000 , 

0.000 

49 

0 

. 146 . 

0.000 

...122 

0.000 

50 

0 

. 17 ? 

. . 0.000 

. 145 . 

0.000 . 

. 51 

_ 0 

.165 

_ , 0.000 ... 

.095 . _ 

0,000 

52 

_ 0 _ 

.195 .... 

0.000 

,113 . _ 

0.000 

53 

0 

...179 

0.000 

.065 

. 0.000 . 

54 

„ 0 

.211 

. 0.000 „ 

_-,. 077 . . ... 

0.000 

55 

0 

.187 

0.000 

.033 . . 

0.000 

56 

0 

.222 

0.000 

.039 

0.000 

57 

0 . _ 

.190 

0.000 

.000 

o.ooo 

5 8 

0. 

.225 

0.0 00 

-.0.00 

.0,000 
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- 

HI- _.1765Er02 

■ Tl-.. 

..6000E+02 

. H2-... 

7i765E-02 

T2-’. 

.6000E + 02 

_ 25 . 

56 

53 

51.. 

. 52 

1 

2 

1 0 

_ -0._ 




' 


. 

EDGE .. 

. 53 

51 


QS- 0. 


-HI- _.,.L7.65E-02 

Jl? 

...6000E+02_ . 

.H2- 

.!i765E-02._„ 

12* 


, 

■■■ 

. 

— _ 

. 







26.. 

_ 56 

.55 .. 

. 53 

56 

1 

2 

3 .0 

0 






— 

EDGE.. 

- 55„ 

. 53. 

.. 

OS- .0, 

— 

_Hl.-„_..m5E.T02 

J !■=.._ 

_tJ6000E+02 H2-_ 

.1765E-02 _ 

T2-_ 

_.6000E+02 

— 27 „ 

58 

57 . 

55 

56 

.. 1 

2 

_1 0 

0. 



. 




EDGE . 

57 . 

. 55 


QS?_0. 


.HI- „.1765E-Xl? 

Jl“™ 

...6000H + 05l1'_ 

J^ 2 

.1765E-02 

T2-. 

__.6000E+02 





1. S .0 ?. ARAMETRIC QUADRILATERAL 

.E.' L_.LM E_J L3 











■ 


.NUMBER OF OUAORILATERIAL eiE«ENTJ_^ 6 

•I::) UMBER OF DIFFSRENT MATERIALS '' « b 


lU3^LQ_e A N e.l_ J.ONJ). 


mater la j. C P N 0 U C T LV I TX 

^liABLE 


J, AY ERS THICKNESS C ON DUCTIVIT Y TENSOR 1 _f HE T A 

too Kjrji: - . 




_.-L0..0PE + 01 .77 1 5_Er.O_3 0.,, ^0 , 

• lOOOE + 01 ^.J715E-03__0. _ 0, 

.lOOOE + 01 ».7 7X5E-r0 3. .0, 0, ' 

.a 0 3 0 E ±0 J ^7 715 £i:i3 J_0 . ^0 . 

.lOOOE + 01 .7715E-03 .0. 0. 

.lOOOE+01 ,1000E-30 0, Of 


ELEMENT .INPUT DATA 


_ . N 



K 

. L .. 






. i 

J 

HATIO 

_.K6.. 

IED6E- 

ISURJ 

. - . 0 






1. 

Wi 

„ 46 

40 

41 . 

. 1 

1.. 

0 

n 

0- 

_ - 2 

A1 

. . 40 

34 

35 

2 

1 

0 

0 

0. 

.3 

35 

34 

2 8 

-29_ 


1 

0 

0 

0. 


29 

28 

- 22 

23 _ 

A 

. 1 

_ 0 . 

. .0.... 

0. 

5 

23 

22 

18 

19 

5 

1 

0 

0 .. . 

0, 

6 

19 

. -.18 

. 14 

15 . 

5 

_-l- 

_ 0 

0 

0. 








— 









... 

— 

— 

— - — 


CO 


o o Jo o o o 



SOLUTION PARAMETER 

S 




TOTAL NUMBER OF EQUATIONS 

= 

58 

SEMI BANDWIDTH. 

= 

11 

NUMBER OF EQUATIONS IN A BLOCK 

Xt 

58 

NUMBER OF BLOCKS 

n 

1 


^TtTK OF 


79 



vjiMi -f' ^ wujrNjf'ji-'y-'i 



f^ODE- NO, N 0 i.A LU-i !i0 +a_V.AiJJ.E N 6 + 2_V.AL.U E_ N 0 ■»• 3 VALU E N 0 +iL _V A L.U Zj 


.308296E + 03 »3 0.2 0 6 A £ +.0 3 . 296.985 E±0 3. . 2938 83.E_+03 , 27.A6 8 lE + 03 

. 2 7„A 2 0 3 E_+ 0.3 25.0J? 03.£ + 0.3 2 5.0m0.£tQ 3 . 2 2 3.8 72.E+ 0 3 223054E + 03' 

.19394 8E + 03 , 19 3 00 3 E + 03 » 180 758 E +0 3 ^,.1739 57E + 03 ..15 82 6 lE + 03 

• ' . 142103E + 03. , 165 5 75 E + 0 3 .1594276 + 03 , 15 4015 E +03 .1429656 + 03 . 

. 15772 3 E + 03. 15.4 08 7 E + 03 . 150.4 2 2E + O3 ^,.1462066 + 03 . 1264 5 3E + 03. _ 

.. .... . .127939E + 03 .152174E + 03 . 149648E+03 .147580E + 03 .144229E + 03. 

• .114 069E + 03 » 115 092 E + 03 .14 9113E + 03 .146887E + 03 . 1450 84 E + 03 

. 1 4206 5 E + 0 3 .103 91.8E.+ 03 ,.104.8 2 9 6+0 3 14 7.4 50E + 03 . 14528 4E + 03 . 

. 14 35 64 E + 03 ..14 0622E.+ 03 . 95 8405 E + 0 2 , 9663 30E + 02_, .14 69 11 E + 0 3 

,1447596+03 , 14 30 53 E + 03 , 1401336 + 03 ,8955716+02 1 . 9 0260 9 E + 02 

.8486836 + 02 .A5.5.054E+.02 ,.8162286 + 02. ,.8221 3.9E.+ 0 2 ,7971646 + 02 

,8028046 + 02 ,7908776 + 0.2 .7.96427E + 02 


03 

o 



i._S_ p.p ,A_R j .R_ X-Q Q U A D R I LA TERAL EL EM ENT X 



CONDUCTION FLUXES 

SURFACE 

FLUXES 

EDGE 

FLUXES 




(LOCAL 

AXES) 

(POSITIVE 

INTO SURFACE) 

( POSITIVE 

INTO EDGE) 



ELEMENT 

OX 

QY 

TOP 

BOTTOM 

I J JK 

KL 

' Ti 












1 

-.2305E+00 

-.3610E+00 

. 0. 

0. 0. 

0. 

0. 

0. 


2 

-.8840E-01 

-.7775E+00 

0. 

0. 0. 

0. 

0. 

0. 


3 

-.3415E-01 

r.fifiS6E+00 

0. 

. 0._ 0, 

0. 

0. 

0. 


A 

-.4254E-01 

-.9997E+00 

0. 

0. 0. 

0. 

0. 

0. 


? 

_-.4353E-01, 

.._r.llllE+0.1. 

0. 

0. 0. 

0. 

0. 

0. 


6 

1913E+00 

-.1194E+01 

0. 

_0, 0, 

0. 

0. 

0, 


7’ 

-.319fiE+00 

-.7861E+00 

0. . 

. 0,. 0. . 


0, 

0. 


8 

.-.2416E+00 

.. -.4072E+00 

0. 

0. 0. 

0. 

0, 

0. 


g 

rf 1522E+00 

-_-.23l2£+00 

0. 

0. 0. 

0- 

0 . 

0 - 


10 

1174E+00 

. -.134BE+00 

0. 

. 0. 

0, 

0. 

0. 


11 

1065E+00. 

7559£rO_l 

0. 

0, 0. 

0. 

0- 

0. 


12 . . 

1066E+00 

-.2464E-01 

0. 

. 0 , 0 . 

0. 

0. 

0. 


' 13 

-.7195E+00 

894iF-01 

. .0. 

. 0. 0. 

0,. 

0.. 

0, 


14 

-..4039E + 00.. 

„-.1085Er0L 

0. 

.0, 0. 

0. 

0. 

"*t?^75E-02_ 


15 

-.2002E+00 

.-.1116Cf00 

. 0. 

-0. ' 0. 

0. 

0. 

-.7344E-02 
-. 7200F-07 

“ 

. „_.16 . 

-.1685E+00 

-.1079E-*00 

0,. 

-0. .0, 

0. 

0. 


17 

1577br00. 

._r.6860E-01_ 

o.._ 

0. 0, 

0. 

n. 



.. 19. .. 

-.1551E+00 

-.23156-01 

■_ 0. 

0. 0. 

n. 

0. 

“-7006E-02 

0 • 


19 

-.5182E-01 

-.7839E <-00 

.0. 

- 0., 0. 

-4 5474E-02 

0 . 

“ 

20 

_.1008E-01_ 

. -.6493E.+00. 

0. 

0. 0. 

-.4639E-02 

0* 

0 , 

“ 

21 

.5616E-02 

-.5251E+00 

0. 

0. 0. 

3qfl7F-o? 

0. 

0 1 


. 22 

.8756E-02 

--4190E+00 

0. 

0. 0. 

. .. -.3461E-02 

0. 

0 • 


- 23 

._.1135E-01 

-.3262E+00. 

0. 

0. 0. 

-43046F-n? 

0. 

0* 


24 „„ 

.1440E-01 

-.244OE+O0 

0. 

0. 0- 

?7?9F-n? 

0 * 

0. _ 

0 • 


25 .. 

.1791E-0J 

-.1696E +00 

0. 

.0.. . . 0, 

R4QQF-ftP 

0. 


26 

_.2199E-0i . 

„-.1006E+00 

0. 

..0. 0. 

-.7350F-n? 

0 * 

a - 


27 

.2678E-01 

-.3489E-01 

0. 

0. 0. 

7?77F-n? 

0* 

n 4. 


- 



. 








o vji ro 


ISOP. A RA METRIC 


0 U A._D .R I. L Jl„T_E..R.. A_L E_X. _E_M ,.E_N ._T_ S 


.CONDUC T ION .F.LUX-.S. SUEFAC^_F_IUX£JS 

- (LOCAL AXES) tP.OSITJ VE—INJO _SUJi£A.C.L) 

element OX QY TOP BOTTOM IJ 


L 06 E_FLUXeS__ 

:eositjve._into.eoge) 

JK . KL LI_ 


1322E+00 

-,9993E-02 0, 

- 0. 

- - 0 , 

0. 

0 

1359E+00 

....... 3011E-01 0. 

0. 

0. 

0- 

n 

1A93E+00 

-,5070£-01 „0. 


_0. 

0. 

_ 0 

2212E+00 

7023E-0L. „ 0, . 

„0 . 

0- 

0. 

0 

3;j02E+00 

1 1A9E + 00 . 0. 

0- 



0. 

0 

01A6E+OO 

-.206BE+00. .0. 

0. _ 

0. 


. . 0 


- - 0 . , .. 

0 . 

. . 0 . 

0 . 




0 . 



OVERALL 


T I H E 


LOG 


N 0 P A_L_ Pp_I NJ_I NPUT s_3_^ 

FORM ELEMENT STIFFNESSES .73 

FORfi TO TAL STIFFNESS ^ 10 . 

li^PO S E BOU.ND ARY CONDITIONS .03 

EQUATION SOLVING 

^LE M ENT F LUXES 5 


TOTAL SOLUTION TINE... 1.47 
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SAMPLE PROBLEM 2 


Nonlinear Analysis of a Convectively Heated, Water Cooled Steel Pipe 


84 





INPUT DATA (SAMPLE PROBLEM 2) 


THFLl 

*«-2 INCH 

PIPE WITH 

water nonlinear 

ANALYS I S 

2/13/1976 


52 

1 

3 

2 

0 0 

1 0 

0. 

0 






1 

1 





0. 


0. 


0.0 


50. 

3 

0 





1 . 


0. 


0.0 



21 

0 





10. 


0. 


0.0 

2 


51 

0 





40. 


0. 


0.0 

2 


* 2 

0 





0. 


1.11 


0.0 


0.0 

22 

0 





10. 


1.11 


0.0 

2 


52 

0 





40. 


1.11 


0.0 

2 

O 

• 

o 

1 0 

25 

1 

0 










1 

0. 

0 

■ 1 










0 • 0584 


0.0 


3 1 , 

0 

2 

0 






1 

1 

3 

4 

2 1 

2 

1 

1.0750 6 

*4940 

10.00 

0.0 

1 .0 


0.0 


0.03472 

300. 








25 

49 

51 

52 

50 1 

2 

1 

1.0750 6 

.4940 

1 0.00 

o 

. 

c 

1 *0 


0,0 


0.03472 

300. 








1 

4 

THERMAL CONDUCTIVITY 

FOR 

TUBE — steel 





32, 


2.208 


212c 

2.167 


572. 

2. 

083 


932 . ■ 

1 .833 

2 

7 

SPECIFIC HEAT 

FOR FLUID-- 

-water 






32. 


1.01 


50. 

1 .0 


1 00, 

0 

.998 


200. 

1 « 00 

300. 


1 .03 


400. 

1 .08 


500. 

1 

.19 




3 

7 

CONVECTION COEFFICIENT FOR TUBE 

TO FLUID 




32. 


0.0561 


50. 

0.0584 


1 00. 

0 • 

064 


200. 

0.069: 

300. 


0.0694 


400. 

0.0670 


500. 

0. 

0614 






PROGRAM OUTPUT (SAMPLO PROBLEM 2) 


_ 25 ..ru f. LL_*12_I N C H_ El P E„W I J H_ W A J EE„ + ♦ _JiDRLI N EAR_AN A UfilL Z./A1/ 13 J.Z 



- - ...... 

1 fL..E_Q._8_M_A_T 1..0 N 





NUHBER.OE_NOOAL POINTS • 5? 


NUnOER.Of. ELEilEMT TYPES • 1 

" ™ - 

— .NUMBER OE..TADLES_ _ - 

‘ - - 

ANAIYSIS..C00E{NANA) - p 

.-ma. .. — 

— EQ.Qj BATA. CHECK ONLY. 


LlnCAK#.. 

'E0,z, ..nonlinear. 

. - 

_PU)T„CQa£(NpL<3T> . - 0 

— — - - 

. - EQ.O. NO PLOTS .GENERATED 

- — ^ ' 

ep.It.unoefdrmeo plot 

- 

EQ.E.._IEMPERATUR£. PLOT 


. _ . MAXIMUM. ITERATIONS - A 

TDI FRANCE » .10000F*00 

^ , 




bTan’k. COMMON LOCATIONS - 109A7 

. - - u - 

.... . _u _ .. — ™ — 

- - ... .... . . 

NODAL POINT ..INPUT DATA 

“ “ * — — • “ 

,, 



NUDE BOUNDARY CONDITION CDDE_ NQDA L .POJ Nt_CQORDlNA TFS 

NUMBER . X V 

Z — KN _ ..TEMPERATURE I 



JL 


„„0 
0 


, 0 -.„ 

-Q— 

_ 0 _ 

0 __ 


Q.OOO __ 

0.000._ 

n.nnn 

0 

<>n.nnn 

_J..000 

O.QOO . 

0.000 

0 

0.000 

.10..QOD 

o.ao.o. 

0.000 

? 

o.nnn 

AQ.aOJ). 

... 0.000 

0.000 

.? 

0# 000 

O-.QOO . 

1.110 

0.000 

0 

0,000 

10.000 

1.110 

0.000 

? < 


-40.000 

1 .110 

-Q.OQO 

2 

0.000 


00 


GENE8ATEtt.N0DAL_QATA 


NDQ£ BQliNBm. 

XDNfi.LLLQN_CHIl£ 

NflllAL_RQIMT COOfiOlMATFS 




NUMBER 

— 

X 

JC 

Z KN 

. _ TEMPERATURE 

1_. 
- ...2 

1 

0 _... 

Q. 

0.000 

0._OQQ..,. 

1,000 

0.000 

L,.U.O 

0.000 

~o7odo 
- Q,ood 

OtOOO 

1 5.0j000 

b.ooq. 1 

0 ■ (too 

4 

5 

6 
7 

0 

D 

0 

0 

1.000 

2.J30.0 

2.000 

3.000 

1.110 
_0,0OQ 

1.110 
0.000 

0.000 

-O.OOQ 

0.000 

o.ooo 

0.000 

.0,000 

0.000 

0.000 




0.000 

^a«.o.ojL_ 

Q..00Q_ 

0.000 
O .Q 0.0_ 

p_.oo'o_ 

o_.po_6_ 

o.boo_ 

_o.ooo_, 

Q.^000_ 

b.6ob_ 

0.000 _ 

o,opo_ 

o.cbo 

0.000 _ 

P..000_ 

o.oop_ 

0»000_. 

0.000„ 

0.000 . 

p.coo., 

p.000._ 

0.000._ 

0,000_ 

0.000,_ 

. 0.000 . 

0.000 .. 

D.000__ 

0.000 , 

0.000._ 

0,000_ 

0,000 _ 

0.000 „ 

0.000_ 

0.000„ 

0,000„ 

0.000 

0.000 . 

0.000 .. 

0,000._ 

0.000 

0.000 

0.000_. 

0.000 

Q..00.0_ 



O.N E_0,L« E.N S..I._O.N_A_L .T_H_E_.R_H_A._L_tJ_L_U„I_0 E_L.E_H_JE..K_I 


^HUMBER OF thermal-fluid. ELEMENTS—' Z5. 

. .NUMBER OF. THERMAL-P.LUID_pftO£E8JJLES_- 1. 


THERMAL-FLUID PROPERTIES [ 

— -T u 8 E_ r7TT^rro3zo ii3zoig"E~]RTf ~ ~ ~ 

PROPERTY CONDUCTIVITY CONOu'cTI V.LT.Y C 0 N V E C T LON CDNVECTION.CONVEC tTo nI__'s7'e'cT fTc'IZ '"spfciEIC VT'^rn’;TTY 

ID K. IA81.E H EJXJP-OMf nr TABLE HEAT. HEAT^TABLE IABLeLI' 


-Qj 1 , 5.8±0 E-0 J CL. 3 . JL DOO £ + OJ 2_.I 1 ^o' 


00 



J)RI(SNAL PAGE IS POOR 


E L C M E fj T. J N P UT DATA 




I. 

a 

'§ 




CONDU CTIO N CO N V,£.C TIO N H A S.S_.F_L OK INLET 


N 

I J K 

L 

PID 

KG 

ISURF 

AREA PERIMETER RATE 

pressure' ' 























1 

1 3 A 

2 

1 

2 

1 

.1075E+01 i65 946 + 01 10_OO.E±0.2. . 

0. 




AF « 

.lOOOE+01 SURFO- 

-O* 


HK- 

.35?2E-01 TFHPK- .3000E+03 HL- 

.3572E-01. 

TEIiP.L 

.JOCOE + 03 


• 







2 

3 5 6 

5 

1 

2 

1 . . 

.1075E+O1 ..65956+01 .10006+02 

.0. " 

- 












— A Fj! 

_ »1000E+01._SU8EQ- 

..0. 


HK- 

• 3572E-01 TEHPK- . 3000F+0.3 HI- 

.3572Er01 

TEHPL? 

. 3000E+Q3 


- 


















3 

5 7 0 

6 

1 

. 2. . 

. . 1. 

-1075E+01 .65956+01 .lOOOE+02 

0. 













AF*' 

•lOOOE+01 SURFO- 

0. 


HK- 

.3572E-01 TEMPK- .30,006+03 HL- 

.. 35.7.2 E-01_ 

TEHPL- 

.3QOOE + 03 




















_ 9 

7 9 10 

8 

1 

2 

1 

.107_5_B + 0_1 .65.95E+01 .lOOOE+02 0. 







* 






AF « 

.lOOOE+01 SURFO- 

0. 


HK- 

_.3A72E-r01 .TENPK- ., 3000.E+03 HL- 

, 3572E-01 

TeHPL~. 

...♦3000E + 03 ._ 









5 

9 11 _ 12 

10 

1 

2_ 

1 

. ... .1075E+01. .6595E+01. .. .lOOOE + 02 

0 . 













AF- 

. ,1000E+01. SURFO- 

0. 


Hlii: 

.3572E-01 TEHPK- .3000E+03 HL- 

.3572E-01 

TENPL- 

.30006+03. . 


















6 

_ 11 13 15 

12.. 

L 

2 

T 

• 1075E+01 .6595E + 01 _. .lOOOE+02 

0 . .. 













AF.» . 

_ .lOOOE+01 SU8F0?- 

o._ 


HK- 

.3572F-01 TFMPK- .30006+01 HI- 

. 3572E-01 

TENP.L- 

_. .30006 + 03 .. _ 

7 

13 15 16 

1.5 . 

' L 

. ... 

1 

.1075E+01 .6595E+01 .IOOQE+0?" 

0. 

AF.^ 

.10006+01 SURFO-' 

0. 



~ hk-..' 

~..'3572E-oT''.TEHPK» "".3000E + 03 HL- 

.35726-01 

T6MPL-.. 

. .30'OOE + 03 ' 


8 15 .17 . 18 16 . 1 , 2 1 .1075E.+01 .6^95F+01 _ _.100J)E.+ 02__0 . _ 

AF« .lOOOE + oi SUkfQ* 0, . . .HK« .3972E-01 TENPK* . . 3000E + 03...HL= . 3572E-01 TEHPL« 7300^E + 03 


o 



9 

_17. 19 

.20 

10 1 . 

2 

1 

_..1Q7,5.E+01_ 6999E+01 ,10003+02. 

0. 

- ~ — w- 

“ - 

— .AF" ™ 

.loooe+oi 

...SORl 






. ..3972E-01 



F0«_.._0. 


HK.* 

, 3 9^72 E=.0 ] ^r.EM£K « 

.30.00E+03 HL- 

TEMPL-.. 

....3000E + 03. . 







— .10 . 

19— 21..„ 

. 22 , 

20 ] 

_2_. 

T 

.L0.15E+01 .6.9.9.9E+01 .lOOOE+02 

_ 0. 














. AF- 

. lOOOE+01. 

_SURFQ« .0.. 


.WK- 

. 3 9 7.2 Er.01 . T.EH P. K • _ 

-_..3000E+03 HL- 

_.3972E-01 

TEMPL- 

* ^non p + n ^ 

' 



... 

— 







_ 11 

21 23 

29 . 

22 1 

_2 

- 1. . 

-1075E+01 ..6999E.+ 01 .lOOOE+02 

0. . . 



_.AF- — 

,.i 9 ooe+oi. 

_SURFQ? 0. 

— 

_HK- 

.3972E-01 .TEH2K- 

»3000E+b3.._HL» 

’ .3972E-01 

TEMPI- 


















.... _ ^ _ 






IZ__ 

23 -.25 .. 

.26 

.29 J 

2 

1 

•.10.75E.+QJL .6999E+01 .lOOOE+02 

0. 










* 



■ 

li 

< 

.lOOOE+01 

SURFQ". 0. 



._..3972E-01._IEMPK-._ 

— .3000E+03.. HL*= 

. .3972E-01 

TEMPL.-— 

— .3000L+03— _ 












— Jia 

25 . 27 . 

28 . 

.26 1 


Jl. 

.1075E+01 .69.99E+01. .lOOOF+0? 

0. 














-AF-„ . 

• lOOOE+01. 

_ SURFQ-. 0. 


JHK? 

..33 7 2 E- 0 ] LEU P.K.- 

. .3000E.+03 HL- 

.3972E-01 

TEMPL« 
























19 

27 ...29 — 

_3 0 

.20 1 

2 

_1 

.1075E+Q1 ,699AF+m ..lOOOF+03 

n. 


1 iHiH 








‘ 



* **" ’ — 

AF » 

_..1000E + 01 SURFO- 0. 


HK- 

•3.9.Z2£=0J LfflEK- 

.30O0E+Q3 HL- 

• 3972E-01. 

tEM.PL- 


























15 

.2.9 31 

.32 

.30 3 

7 

1 

lL0l 7_5E.+ O1 _,6999£ + 01 .lOOOF+rs? 

0. 













' 

AF? _ 

.lodoEtoi 

SURFO- 0, 


HK- 

,-3972E.-01_ .TEMPK- 

.30OOE+03 HL- 

..,39726-01' 
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MATE.Rl. AL PR Q .P_E..R.„T._J_£_.S .T_A 


J a'blE_NUMBE R £ ] JjH tR HA UXQMLUj;XlVJXt_FJlR_TmE ™ itEEJ. 


XE M P • MD E E RXY IE.tt.P..« ER.OEE .RTY TEM E* PRo'p'eRtV T EMP." 'VR nVpR'TV 


32.0 

.2208E+01# 

« 

.212. 0__ 

_.216.7E + 01> 

5.72,0. ..2083E+01, 

,_932.0 

. 1833E + 01; 







1 








table_ 

.NUMBER. 2 

SXX.C.I FI CHEAT F □ R 

FLUTD — WATER 











_ X£MP._ 

P.SOE£g_IY 

TFMP- 

PROPERTY 

TEMP. PROPERTY 

TEMP. 

PR nppR rv 


32.0._ 

.. .1010E+01> 

—50...0 . 

lOOOE+01 . 

100.0 .Q<38nF+nn, 

200.0 

. lOCOE+Oi r 

3QQ.Q 

_,1.03 0E+01> 

.APJL. 0 

_.,....1.0 8 0E+01^ 

500...0 .1190E+01, 



. 


1 

table„ 

.NUMBER 3, 

C aiM.VXCJJ D H_C.D E FXI C I.E N.I_JE£lR._r U B E_.T 0„.Fi UU3_ 




„ ._J E M P . PROPER Txl J.Eiii. ER.O P.EBxVI IE ROE ERXY_ I.E M.E.. PROP ERJ.y. 


32.0 .5610E-01> 

50.0 

,5.8A<2£r:.0O 

100.0 

.6A00E-01 , 

200.0 .6930F-O1 

.300. 0_ .69^0Er01, 

iO.0. .iJ 

_..6XO.O.E-OLje 

5.0.0 ».0 

,.61^0E-0JLj* 
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SOLUTION PARAMETFR 

S 




TOTAL NUMBER OF EQUATIONS 

- 

52 

SEMI BANDWIDTH 

s 

4 

NUMBER OF EQUATIONS IN A BLOCK 

- 

52 

NUMBER OF BLOCKS 

s 

1 



94 



X ..E_H _ P__E _ R _A_L„U_ R_E 


V ECTOR 


_NON L £AJR_A N ALYSIS ITFRATIO N NUHBFR 1, 


6 P L.N Q, NO ]/ 1 L U.’t N[0Jil„V.ix.aE N G + Z VA LU E N O + 3 ~VAUJ. E.__... ", N D + A VA L u'e.7 


1 .50000 0'= + 02 .15.A050E+03 ,33.3 261E + 0Z ,154A57E + 0 3 .5665 3^E + 02„ 

.6 , 155^50E + 03 « 59.a.2 51E±02 ..15.67_B8E+0X ^._630282E+.0? , lt>8325E + 03.__ 

'll . _ . .660914E + 02 .159973E+03 . 691989E+02 -161680E + 03 .721716E + 02 „ 

-.16 . . . .163912E + 03 .751953Ef02 .165151E+03 .780855E+02 .166885E + 03 

21 .810 310E + 02 ^,.16.8 6.08 E + a3 ,.8.668 3 5.E±02 « 172001 E + 0 3 .922 A 65 £ + 02 

26 . 175318E + 03 ^976060EJi02 . 178 5 99 E +03 ..1028 87E + 03 .181700E + 03 . 

31 .107969E + 03 .189767E + 03 ..112982E + 03 .187757E + 03 .117801E + 03 _ 

36 . 1906 6 7E + 0 3 ..12 2 5 60 E + 03. U 9 3 5.0 3 E .+ 0 3 12.7X2 9E + 0 3 , 1962 60E + 0 3 

91 . 131697E + 03 .198992E + 03 . 135977E + 03 .201539E + 03 , 190263E + 03,.. 

^96. . 209019E + 0 3 . 19 9 35 7E+0 3 . 206296E + 03 . 19 8 397E + 0 3 .208172E + 03 

5.1 . 1 5 2 1 9 8 E + Q 3 ^ZO.3.0^1 £ + 0 3 


- LARGEST.CHANGE— ..LO.O..QO.O 0 E+ 0.3-.&£ R_C E N.I_ Ll.Ji 0 P.E Z 


^o 



TEHPERATU R E 


VECTOR 


Q 




ss 
§8 

ss 

|E 

tHi a 

§ 

o 

O t-3 
£d a 
a 


NONLINEAR ANALYSIS ITERATION NUMBER 


NOD E NO. 


NO 


VALUE 


NO+1 VALUE 


NO+? VALUE 


N O-t-3 VALUE 


_. 1 \ 

11 
16 
. 21 . 
-2 6 . 
J 1 
.3 6 ^ 
41 
46 . 


_.500000E + 02_ 
..,14a5 5 6E_+.0_3_, 
.672946E+02 
.156340E+03 
, 634842E 'f02 


.. 147_18 0E±03. 
-♦-60..5.3 7XE+_02_ 


-U4.9_8A0X-ta3_ 


.. .152968E + 03 „ 
„ ♦771548E+02 
_*_161A8 6E±iI3__ 


706636E+02. 
_ . 158056E + 03 
_._ 8 . 9 i >_ La 6 -E + 02 . 


168 2 9.4 E±.Qi_ 
_.112867EtQ3_ 
e_49l7E + 03_ 


-♦_63.9.999E + 02- 
..154639E + 03 
.,802763E+02. 
_..1.64912E.+ 0 3_ 


,138451Et03 - 

X.. .199918E + 03. 

^ a 6044 , XE ± S .3 


.jJ.Q.15_Q.9,E.+53_ 
_!l.178 3 5.3E + 03_. 
.jJ.2.8_6.4.0_E+.03_ 
194190E + 03, 
_.152095E.+ 03_ 
_^.0_lZI2.a±Q3- 


...17.16.3JE+_0.3_ 
...118_2_61E + 03. 
35.08.9Et.03_ 
143103E + 03, 
_».2O2 5 30Et.03. 


.^_l-0.Z2,.8.1E+.0 3_ 
_.i.l81673E + 0 3. 
_^_1335 6.9E.+0.3_ 
_a97109E + 03. 
„...156424E + 03_ 


largest CHANGE._.5.1_41430E+01 .PERC ENT._A IL.MDDE__„ .5 1 


NO-(-4 V A1.UF> 


.,571364E + 02 
-.151359E + 03 
.738686E+02 
-,159773E + 03 
_*_956770e + 02 
_,174967£ + 03 
...123480E + 03 
_,19U78E + 03 , 
.♦147717E + 03 
.•204720E1 03 


to 

O' 




. R. A T U R E V E 

-C- T Q R 







NONLINEAR ANALYSIS 

TTFRATION NUMBER 

3 

** 




NOD£_bj.O..__ 

- NO VALUE 

tlQ+1 VALUE 

NO+2 VALUE 

NO+3 VALUE 

Nn+4 VAI UF . 



. 1 . 

_ „ .500000£ + 02„ 

.lA6779E+03_ 

._535620E + 0.2._ 

.1471 84F+03 

• 571547E + 02 

-6 

..♦148157E + 0.3__ 

..6.CL5 6.A2E+.02 

Ui!9A6.2E+03_ 

.840361 F+0? 

♦ L50963F+03 

11 _ . 

_ .673390E+02 

.152573E + 03 

.,707179E + 02 

. 1 <54?48F + 03 

. 739310E + 02 

_ 16 . 

« 1559A9E+03 

.77Z271E+02 

. 157667E+03 

. 803869F+0? 

, 1 5938 8E + 03 

21 

..8367A6E+P2 

.161106E+03 

897276E + 02 

. 1 64580F-5-03 

. 968045F + 02 

._.2 6. 

.167969E+03 

. 1 01 858E + 03 

. 171395E+03 

. 1 074ZqF+03 

.174838E+03 

31 _ . 

- .112958E+03 

.178238E+03 

. 118416E+03 

.181 578F+03 

. 123638E + 03 

36 

.18‘i83ZE + 03 

.128802E+03 

« 188Ql^E+03 

.'1 33734E + 03 

,191I13E+03 

‘^1_ _ 

, 138618E^-03 

_..19A135E + 03 

. 1A3272E+03 

, 197053E+03 

• 147889E+'0 3 

4 6 

■ .1998B1F+03 

.15226BE+03 

. 20250.2E + 03 

. 156598E+03 

, 204699E + 03 

51 

♦160621E+03 

.-2IL525J5£+D.3 
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LARGEST CHANGE._„ .2.I29551E + OQ P£R.CEhiT.„AT._NQDE 
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TEMP ERA T U R E 


V E C T □ R 



NONLINEAR ANALYSIS 

ITERATION NUMBER 

4 



» 


NODE NO. 

NO VALUE 

NO+1 VALUE 

NO+2 VALUE 

NO+3 VALUE 

NO+4 VALUE. 


1 

,500000E.t02_ 

. 146758 E + 03 _. . 

..5 3 5623 E.±.O.a.. 

- 197163F+03 

.571557F'*-02 

6 

.148136E+03 

.605657E+02 

.J49A4 2£jL03 

.640380F+0? 

. 1 509A3F-f 03 

11 

. . .673414E + 02 . 

... *152553E + 03. 

... 707207E + 02 

.. . ..154226E + 03 

. 739343E + 02 

16 

.155930E+03 

.772309F+02 

157649E + 03.. 

•803610F+02 

.159370Ef03 

21 

.835792F+0? 

. 161089E+03 

. 897330E+02 

. 1 6A535F+03 

. 958107E + 02 

26 

. 167957F + 03 

.101665E+03 

.„171388£ + D3 

- 1 n7A36F+03 

. 1 7A835F+03 

31 

.112965E+03 

.178237F+03 

..118422E + 03 

. 1 81 575F+03 

. 1 236A5F + 03 

36 

.184832E+03 

.128809E+03 

.1880I5E+03 

.1337A0E+03 

. 191113E + 03 

.41 . . 

_ _.,138625E + 03 

. 1.941 36F + 03 

. 143279F + 03 

. 1 9706AF+03 

. 147895E + 03 

46 . . 

.199683E+03 

. 152274E+03 

.202.503E+03 

. 15660AF+03 

« 204701E + 03 

51 

. 160627E+03 

-205757E+03 
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1.4 4 0 7 6.9 E -zOX. P E R C E N T_A T_KO O.E 
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C DN V E R G E NC E AC HIEVED 


00 



„0 H i." 0 *I H^E. « rS_r_OJl_A_l 1_U_E_R._H_^L_-^_L_U.JLJ3 e-JL_e_B_E_N_T’'s_’’_.' J J _ZIIIZIZTIZZZZI^^ 


„.H. e A_T. E„LU„X_£-S f_L.U. I_0_. P. R.E S_S_U.R. £ ._O.X-L A„, 


EL6HENT .. 

. FLUID . 

TU9E 

EXTERNAL 

FRICJION P.RESEURE 

now ACCELERATION. 

. .PRESSURE 

. PRESSURE 

HEAT FLUX 

.CONDUCTION 

. CDNVECT.ION FLUX 

DROP,, . . 

PRESSURE DROP,. 

__ NODE-I.. 

NOD£-J._ . 














1 

.5170E+03 

-.9A91E+00 

.3951£tOZ . ... 

n i. 

0. .. „ 

.,0... . 

.0. 

2 

.5535E+03 

. -.2283E+01 

.3A35Ef02 

0. 


0. 

0. 

3 

.5eaAE')^03 

-.3Q61E+01 

.. .3A09Ef02 „ 


„ , 0. 


n. 


•6227E+03 

-.3519E+01 

•3378E+02 

. 

0. . _ . 

0 . . 



.6565E+03 

_-t3775E+01. 

• 33A3E+02 

0.. 

.04 

.. 0 . 

0. 

- 6 

.6896£t03 

-.3922£*01 

. . 3306F+02 

0- 

, n- 

0. 


7 

.7226E+03 

-.3992E401 

. .3265E+02 


0# _ , . 


- -- 

^ _ 0- 


•7551E+03 

=.9Q29£+01 

.3229E+02 

0. 

0. . 

0. 


9 , 

•7371E*03 

“.A031E+01 

...,3190E*02.. . 

.0.,...,. 


„ 0. 

.0.. _ . __ 

10 

.ei87E*03 

_-.A026EF01 

.3151Et02 

0. 

0. 

0. 

0- 

11 

..8653E+03- 

-.A03‘iEt01. 

•6186E+02 

0. 

0. 

0. 

0. 

12 - 

,9261£*03 

_-.9005£+01 

•6032F+02 

0. 

0- 


0. 

13 . 

♦9855E+03 

-.AOlAC+01 

_ .5877£*02 

0.... .. 

0. 

0. 

0. 

} A 

.lOAAEtOA 

„ -.'.031E+OL 

.5722Ct02 

0. 

0. 

0. 

0, 

15 

. IIOOE+OA 

. -.3977E«^01 

.556flE+02, 

0. 


0. 

n* 

16 . 

.1155E+0A 

-.3901E+01 

_.5A 16E+02 .. 

0. 

rt . 

n* 

n* 

17. 

..I2oee+OA. 

.-.3805E101 

,5257E*02. 

0. 

0. 

0. 

0. 

. 15 .. 

.1260E+OA 

.-.3718E+01, 

.5122E+02 

0. 



0. 

-19 

.1311E+0A 

.-.3618E+01 

.A93nF »n? 

. 

0. 

0. 

0. 


_.1360E»0't. 

-..352BE + 01. 

.A8A2Ef02 

f) 

0. 

0. 

0. 

21 

.. IAOOC*OA. 

=.3A16E+01 

.A708EF02 

0 . 

0. 

0. 

0. 

_ 22 

. 1A5AE+0A 

.-.3288E+01 

.A578E+02 

0. 

0. 

0. 

-0. 

23 

..lA99Et09. 

_ -.3056E + 01. 

.AA56E4-02 

0. 

0. 

0. 

0- 

. 2A 

•15A3E+0A 

-:.2562£ + Ql 

... .A3A7E4-02 

n. 

0. 

0- 

n- 

25 

.1585E+0A 

-.1231E+01 

• A27AEF02 

0. 

0. 

0. 

0. 



OVERALL TIME LOG 


NODAL POINT INPUT... 

.15 

FORM ELEMENT STIFFNESSES. 

. A9 

FORM TOTAL STIFFNESS 

.93 

IMPOSE BOUNDARY CONDITIONS 

.06 

EQUATION SOLVING 

.lA 

ELEMENT FLUXES 

.15 


TOTAL SOLUTION TIME 

1.93 




eeproducibility of JPB 

TOR 
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SAMPLE PROBLEM 3 


Nonlinear Analysis of a Simplified Heat Exchanger 
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INPUT DATA {SAMPLE PROBLEM 3) 


SlMPLlFltO HEAT EXCHANGER ■X#EAT*'«' K*NQNL 1 NEAR*#MEMOW I AL DAY. 1976 
51 1 6 2 0 0 


1 

1 

0.0 

0,0 

0,0 

3 

250, 

49 

1 

40.0 

0.0 

0.0 

3 

250. 

2 

I 

0.0 

0.75 

0,0 

3 

50, 

5 

0 

2.5 

0,75 

0,0 

3 

0, 

50 

O 

40.0 

0.7S 

0.0 

3 

0, 

3 

1 

0.0 

1 .50 

0.0 

3 

500, 

51 

1 

40.0 

1 .50 

0.0 

3 

500 » 


11 16 1 0 INTERNAL FIN 

t 0.0231 6. 


0.2513 0.0 

1 

0.7 

2 3 





114 5 

6 

3 2 

1 3 

0 

41,04 



0,3 1.0 

1 ,o 







16 4 6 49 50 

51 

48 47 

1 3 

0 

41,04 



0.3 I.O 

1 , 0 







1 1 2 FLU 1 D 

TO WALL 

convection 

COEFFICIENT 

~H tBTU/FT2SF ) 


40.0 0.2513 

46.80 

0.2627 

53,53 


0*2737 

60. la 

0.2848 

66,77 0,2956 

73,29 

0 ,3066 

79,75 


0*3175 

86*15 

0.3284 

92.49 0.3397 

98.77 

0.3S01 

105,00 


0*4011 

105. 

0.401 1 

2 6 SPECIFIC HEAT 

“ C SUB P (BTU/LBM-Fl 





“50.0 0,63 

0.0 

0,675 

60.0 


0,735 

120.0 

0*B05 

140,0 0,84 

220.0 

0,850 






3 20 FLUID 

VISCOSITY - MU (LSM/FT“SEC) 





"-50. 0,2450 

“4 0* 

0, 1 350 

“20,0 

0 

.052 1 3 

0*0 

0,02352' 

20. 0,01220 

4 0, 

0 .007299 

60,0 

0 

,004649 

80*0 

0.003049 

100, 0,002081 

1 20. 

0.001581 

1 30*0 

0 

,001395 

138*0 

0*001269 

140. 0,001243 

142* 

0*001219 

160,0 

0 

.0009878 

180,0 

0,0007849 

200. 0,0006248 

220. 

0,0005131 

240,0 

0 

,0004394 

300*0 

0*0002960 


A 11 FRICTION FACTOR F 

40.0 0,01220 46.80 0.0II92 53*53 . 0-01161 60*18 0.01133 

66.77 0.01109 73.29 0.01174 79,25 0*01142 86-15 0.01109 

92,49 0.01072 9B.77 0,01040 105*00 0-01012 

5 O FLUID density - RHO (LBM/FT3 ) 

-50. O ■ 70.0 0.0 69-1 50,0 60-2 100, 67.15 

150.0 65*8 175,0 65,0 200.0 64,25 225. 63.5 

6 4 thermal conductivity, K- {BTU/SEC-FT-F 1 

32.0 0,02194 212. 0,02306 392, 0.02444 


o 

04 


572 


0*02503 
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PROGRAM OUTPUT (SAMPLE PROBLEM 3} 


1 IiLP.LLEJ.£b„H£ A I^dgjaAUGEA ♦♦IaJ** **RQMLXM£A&1*H£Mq\i ALD Ay!>1'i.Q76 


• .C._D„N_T_R_D_I I NFO R H A T I n N 


-NUHBER;_0f._N00AL POINTS ■ SI 

NUMB E£_0£_e I, eH E NJ_.TJt£.E.S - 1 

NUMBER ,OF„TABLES 6 

ANALYSIS. CODE (NANA) 2 

E Q .Aj> OAIA_CH E£J(_0 NJJtj 

:_E(3.1i .LINEAR^ 

E0.2* NONLINEAR ... 

P.LOT_CDO£(NPLQn - 0 

.. EQ.O# NO PLOTS GENERATED 

EO.l/ UNDEFORMEO . P.LOT. 

E Q , 2jl_.LE M P.ER AIUB E_a 0J[ ^ 

I.TERA.TION_EARA«E.TERS ! 

MAXI HUM _I.TERAJ.IONS - 6_1 

TOLERANCE - .1 0000E4.Q.Q 


RJ.AN K_C.QH MON_L OC A.T.I OH5_s_10l9.6.J 


NODAL 

POIN.T INP.UT..QAJ.A 













J(DDE — 

BOUNDARY CONDITION CODF 

NDRAL-LOLNI 

.CDIIfi.DiHAr£S_ 




.NUMBER 


- X 

Y 

Z KN 

TEMPERATURE 








1 




D.DDQ 

■K 

?5n .non 

.A.9„ 

1 

AO.OOO 

WtSKSIStttK^M 

Q^QOQ 

3 

250.000 

a_ 

1 

0.000 

.750 

0.000 

3 

50.000 

5 

a ' 

2. boo 

* 7*50 

0.000 

3 

1 0.000 

50,_ 

Q 

AO. 000 

.750 

(LiiLOO 

3 

0.000 

3 

L_ 

0.000 

1.500 

0.000 

3 

500.000 

5J_ 


40.^.000 

-J.SiLO 

_a.^OJDCI 3__ 

500.000 


— — - 

_£EN£RA.T ED -NODAL DATA 



■ nil Miiim iiiiiiiM 


NODAL POINT 




JIUMMR- 



X 


I 

KN TEMPFRATURF 


1 

1 


0.000 

0.000 

0.000 

250.000 

? 

1 


0.000 

.750 

a. ooo 

50.000 

3 

1 


0.000 

1.500 

0.000 

500.000 



1 





750.000 

_5 

a 


2. .50.0. 

.750 

O.OOQ. 

0.000 

. 6 

1 


2.500 

1.500 

0.000 

500.000 

7 

1 


5,000 

0.000 

0.000 

750.000 








o 

cn 


e 

0 

5,000 

.750 

0.000 

0.000 

9 

1 

5.000 

1.500 

0.000 

500.000 

10 

1 

7.500 

0.000 

0.000 

250.030 

11 

0 

7.500 

.750 

0.000 

0.000 

12 

1 

7.500 

1 .500 

0.000 

500.000 

13 

1 

10.000 

o.coo 

0.000 

250.000 


l<t 

0 

10.000 

.750 

0.000 

o.oco 

15 

1 

10.000 

1.500 

0.000 

500.000 

16 

1 

12.5 00 

0.000 

0.000 

- 250.000 _ 

17 

0 

12.500 

. 750 

0.000 

0.000 

I'a 

1 

12.5 00 

1 . 500 -_ 

0.000 

500.000 

19 

1 

15.000 _ 

0.000 

. 0.000 „ 

250.000 


0 

15.000 

.750 

0.000 

0.000 

21 

1 

15 , 000 .. 

_ _ -. 1.500 

0.000 

. 500.000 

22 

1 

17.500 

0 , 000 . 

0.000 

_ 250.000 . 

23 

0 

17.500 

.750 

0.000 

0.000 

2 ^ 

1 

17,500 . 

1,500 

0.000 

500 . 000 . 

25 

1 

20.000 

0.000 

0.000 

250.000 

26 

0 

20.000 

.750 

0.000 

0.000 

27 

1 

20.000 

1.500 __ 

0,000 . . 

500.000 

28 

1 

22.500 

0.000 

0,000 

250.000 .. 

29 

0 

22.500 

. 750 _ 

0.000 

0.000 

^0 

1 

22.500 

1.500 

. . 0 . 000 . 

500.000 

31 

1 

25.000 __ 

J 0 . 000 .. 

0.000 

250.000 


32 

0 

25.000 

,750 

0.000 

0. 000 


33 

1 

25.000 

1.500 

0.000 

500.000 . 


3 ^ 

1 

27.500 

0.000 

0.000 

_ 250.000 


35 

0 

27.500 

.750 

0.000 

0.000 


36 

■) 

27.500 

1.500 

. 0 , 000 . . . 

500,000 


37 .. 

. .. . 1 . _ . 

30.000 

0 . 000 „ „ 

. 0,000 

250.000 

2 ^ 

38 

0 

30.000 

.-750 

0.000 

0,000 

39 


30.000 

1 , 500 . 

0,000 

500.000 

Ss 

<vo 

1 

32.500 

0.000 . 

. 0 . 000 . 

250.000 

^ 5 


0 

32 . 500 . 

- ,7 50 . 

0.000 

0.000 _. . 

is 

92 

1 

32.500 

1.500 

0.000 

.. 500.000 .. 

tel B 

93 . 

1 

35.000 

0.000 

0.000 

. . 250.000 

99 

0 

. 35,000 

.750 

0 . 000 . 

0.000 

95 . 

1 

35.000 .. 

-_ 1.500 

0.000 

.... 500.000 

96 


37, 500 _ 

.. 0,000 

. . _ 0.000 

. 250.000 

47 

0 

37.500 

-.750 

0.000 

. 0,000 __ 

_ 98 

1 

37.500 .. 

..„ 1 .500 

0,000 ..^ 

500.000 _ 


4 9 

1 

40.000 

0,000 

. 0.000 , J _. 

250.000 

50 

0 

40.000 

.750 . . 

0.000 

0.000 

o 

w w 

tel 

.51 

1 

40.000 

1,5 00 

0.000 

. 500.000 
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TWO.OIMEH S_1_0 M A. L T_H..E_R. M A _L-_ 0 E..L_E H_E.A!.T . _ 


_NUMBER OF THERMAL-FLUID ELEMENTS.- 16 

_HUMB6R OF THCRHAL-F.LUI0_P_R0PER.LL£3_» 1 


.„INtEJ^NAL„EIN. 


THE R M A L - fI U I D~ P RO P E R T I~E S 


- N E_L.., U ...L-D R_jR_0_£-J_S_t_JLJ-_S. 


PROPERTY CONDUCTIVITY CONDUCTIVITY CON VEC TI01i„ CONVECTION .CONVECTION SPECIFIC SPECIFIC VISCOSITY.. 

.. _J0 K T.AB.LE H EiEOHENJ LABLE KEAT HEAT. .TABLE TABLE 


._.23i0E.-0.X- 


.j.^L3E±0iL 




_^O.OO.E+O.0„ 


3 


ss. 


S3 



EUMENT 'im>0t'6ata 


FJN FID area TOP ^BOTTOM KASS_FLOW INLET 

. £ 1 - I- K I H N PIP K G THIC K NESS E FFICIENCY F^Q.TOR. Wl DTt! WIQIH RAJS PP’Is SUP e!! 


1 

1 


5 

6 

3 

2 

1 

3 

.3000E+00 

.lOOOE+01 

lOOOE+01 

.7000E+00 

.7oodE+do 

.AlOAE+02 

0* 

2 

A 

7 

8 

9 

6 

5 

1 

3 

•3000E+00 

•lOOOE+01 

. lOOOE+01 

•7000t+00 

.70006+00 

10AF402^ 

0, 

3 

7 

10. 

„1 1_ 

.12. 

9 

8 . 

.1 

3. 

.30O0E+OO 

.lOOOE+01 

. lOOOE+01 

•7000t+00 

-7000F*fjO 

. 4l O^F+0? 

0 , 


10 

13 

... 10 

15 

12 

_..J.l 

.1- 

.3 

.3000F+00 

.1000E + 01.._ 

.lOOOE+01 

, .7000F+00 

70006 +00_..’' 

.AlOAE+02 

0.' 

5. 

13 

JO. 

. .17.. 

18. 

.15. 

lA 

.1. 

3 

.3000E+00 

.lOOOE+03. 

.. lOOOE+Ol 

. .7000F+00 

.7non£+oo 

•A1CAF40? 

n 

6 

16 

19 

20 

21 

18 

17 

1 

. 3 _ 

. . .3OO0EfOO . 

•lOOOE+01 

.loooe+oi 

.7000E+00 

. .70006+00 ” 

.AlOAE+02 

0.' 

7 

19 

22 

23 

2A 

21 

20 

1 

. 3 

.3000E+00 

.lOOOE+01 _ 

. looor+oi 

.70006+00 

. 7000F+00 

.AIOAc +02 

0. 

0. 

22 ... 

25 

.20. 

. 27 


23 

J 

3 

.. 3000E+00-. 

..JOOOE + Ol 

. looor+oi, 

. _.7000C+00 

. 7000E +00 

• '•10^^40? 

0 , 

9 

.„25 „ 

23 

. 29.. 

™ 30 

_ 27 

26 

1 

3 

. 3000E + 00, 

.lOOOFtOl 

. lOOOL+Ol 

_ .7000C+00 

_. .7CCCE+00 _ „ 

.AlOAF + 02.' 

"o." 

__10 

26 „. 

31 

32 . 

_ 33. 

30. 

29 

1_ 

3 

.3000E + 00 

,i000t+01 _ 

.lOOOF+Ol 

.7000E+00 

„.7000fc + 00 _ 

♦AlOAE+O? 

0.. 

n 

. _J1 

3A. 

35_ 

36 33 32 

1 

3 

jJOOOE.+ OO 

.t.lOOOE + 01. 

.. lOOOE+01. 

. 7000C + 00 

• .70001. fOO 


0* 

. 12 

3 A 

37 

. 38 

. 39 

36. 

...35 ... 

1 .. 

3 

_,3000E+00 

.lOOOE +01 

. lOOOE+Ol 

,70006+00 

_. .70006+00 .. . 

.AlOAE+02 

0. 

13 

37 

AO 

. A1 

A2 

39 

.. 38 

l.„ 

3 

.3000E + 00 

.lOOOE+01 _ 

.lOOOE+Ol 

.70006+00 

. 70006+00. 

.A10A£+02_ 

0. 

- i ^ 

- AO 

03. 

AA_ 

__A5.. 

.02. 

iJ, 

,} 

3 

,3000£+O0 

,1000E+01.„ 

.lOOOE+Ol 

.70006 + 00 

70001 +00 

.AlOAE+02 . 

0. . 

I’j 

A3 - 

A6 

A7 

. A8 

. A5 

-_AA .. 

1 _ 

._3 

, 3000E+00 

.lOOOE+01 _ 

.lOOOL+Ol 

.70006+00 

.7000E+00 

.AlO+E+02 

0. 

10 

A6 

A9 

50 

51 

AB 

_ A7 

1 _. 

_.3 . 

.3000E+00 _ . 

.lOOOE+Ol _ 

.lOOOE+01 

.70006+00 

..7000E+CO . 

.AIOAl+02 

0. 
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H A T E R. I A .L „ . P R. n P F R T r P 


TABLE NUMBER 1 FLUID TO WALL CONVECT ION COEFFICIENT -H( 6TU/„F„T2SF L 


T 6«P. PROPER T.Y T E M P_» £ R.Oiei T Y llE PR 0 P_E RT Y~ t>mV. pg q P,E P T Y 


AO.O 
66.3. 
92. f) 


.2^>13E^OO^ A6.8 .2627E+00# -53.5. .2737E+00^ ..60.2 , .2BA8E--00, 

,2956FtOO/ 73.3. ,3066FvOC)j' 79,3 .3175E+00;. 86.2..._ ,3^8AE^00/ 

.3397E+00# 93.8 .... 3501E+00/ J05,0_. . AOl IE +00.». 105.0 .AOllE + 00/ 


. . „TABLE .NUMBER ...2 SPECIFIC HEA T - C SUB P (RT1J/| BM-Fj 


lEMP. . PROPEP,TY TEMP.. PROPERTY TEMP^ P.ROPE.RTY .TEMP.. T'pROPERT'y 


-50.0 . .630CE+00.P 0,0 .,6750E+0Q> 6.0...0 ,_I350E.+QQ> 120. 0__ .eOSOE + OO# 

lAO.O . .8600E+00^ 220,0 .8500E+00.. 


. . .table NUMPE.R 3 LlUJD VISCns'lT Y'~-"~M U (L BM/FT-SFC ) 


T-E.M P_., p P 0 P.EAT Y ZlEM.P . PkaPE RTY TEMP. PRO P ERTY T EMP. PR.p P.E R.f t. 


— -50,0 .._..2A50E+00,_ -AO.O ^..1 35OE + 00.T r:2 0 . 0 .5213E_-.01> 0.0 .2352E~-01^ 

PO.O . 1220E-01> A 0.0 i.7299E-02» 6 0.0 , A6A9E-02^ 80.0 .1, .30A9£-02^ 

- .100.0 .2001E-02. . . 120.0 _..,1581E-02, 1 30 , 0 . 1395E-02V1 ik 8 . o'l _Vl2 69E-02/ 

190,0. .12A3E-02/ 192.0 _.1219E-02# 160.0 ,9878E-03» 180.0 ,7899E-03^ 

200.0 .6298Er03/ 2 20.0 .5.1 31 £-03;. ,29.0.0 .9399E-03^ 300.0 .2960E-03/ 


TABLE .NUMBER „ 9 , .F.RICTION .FACTOR _- £ 


- _TEMP. _ PROPERTY TEMP* £.RPPE.R.IY T E'm'p . PR.Q£E£IY. 


JEMP. .PROPERTY 


-.{'P-.-Q •.! 22.?Er.0_li ..I 192 E-Q1? 53.5 60.2 . 1 133E-01^ 

_ . 66. U..., ,li09E-01j._ .73.3 ^.],179E-01^ 79.3 . 1 1926-01/ _ ~q’672 .’u’o9r~01/ 

92,5 . 10726-01 . “(i.P .1090E-01/ 105.0 .lOlZF-Ol/ 



8LE 

NUMBER 5 

FLUID 

DENSITY - RHO 

(LBM/FT3 ) 

/ 












TEMP. 

PROPERTY 

TEMP. 

PROPERTY 

TEMP. 

PROPERTY 

TEMP. 

PROPERTY 


-50.0 

__jl700PE + 0 2j! 

0.0 

.6910E+02. 

50.0 

.6820E.+ 02, 

100.0 

.6715E+02, 

... 150,0 _ 

65806 + 02 ^ 

175.0 

.6500F+02, 

200.0 

.6A25F+02, 

225.0 

♦6350E+02, 









• 









TABLE 

NUMBER 6 

THERMAL CONDUCTIVITY 

, K-(BTU/SEC-FT-F ) 











IEMP.._ 

PROPERTY 

-jjEMpr 

PjjoVLRTY. 

"lemp. 

PROPERTY 

TEMP. 

PROPERTY 


32.0 

__.21.9AE.-0_1, 

212.._q_ 

.,_2.3.0.6E-01.,. 

3.9L2...D 

..2_AA„AEi^ai, 

5JZa.Q 

.2583E-01, 



a 


mT_A.L_mJJl BER OF EQUATIONS 
S£ M I_ B A N Q. WJLQXH 


5JL 

-b. 


NUfiBER QF..EQUA' 


N U M B E R OF R i nCKS 


1 


OF 

CfT055^ _ ^r\ T>nOiB> 



TEMPERATURE VECTOR 











NONLINEAR ANALYSIS ITERATION NUMBER 

1 



NODE NO. 

NO VALUE 

NO+1 VALUE 

NO+2 VALUE 

NQ+3 VALUE 

NO+A VALUE, 


1 

.250000E+03 

.500000E+02 

. 500000E+03 

.2500006+03 

.6030296+02 

6 

.500000E+03 

.250000E+03 

. 702241E+02 

.5000006+03 

.2500006+03 

. 11 

.79B882E+02 

. .500000E+03 _ 

.250000E+03 

,„..89i895E + 02. .. 

,5000006 + 03 

16 

.250000E+03 

.9825^7E+02 

. 5000006 + 03 

.2500006+03 

.106975E+03 

21 

. 500000E+03 

.250000E+03 

. 1] 5478E + 03 

. 500000E+03 

.250000E+03 

26 

. 123653E+03 

.500000E+03 

. 250000E+03 

. 1 316306 + 03 

,5000006+03 

31 

.2'30000£ + 03 

.13929AE+03 

.5000006+03 

, Z50000E + 03 

.1A6777E+03 

36 

. 500000E+03 

.250000E+03 

. 153961 F+03 

.5000006+03 

.2500006+03 

A1 

■ .160981E+03 

.500000E+03 

. 2500006+03 

. 1677156 + 03 

.5000006+03 

A6 

.250000E+03 

.17A301E+03 

. 5000006 + 03 

.2500006+03 

, 180613E + 03 

51 

. 500000E + 03 



















LARGEST CH AN6£._..,1 000000 E + 03.. P_ERC ENT.. AT...NODE 5 

I 




.637866E+02 
♦ 2!50000E + 63 . 
.500000E+03 
, 13235^6 + 03 , 
,2500006 + 03 I 
,5000006+03 ’ 
. 187^166 + 03 
,250000E+03 ' 
,500000E^ 03 
.2298^16+03 


LARGEST CHANGE__..2176A5,8£+0Z .P_ERCEKT._AL_N00.6 



1- . ,2500006 + 03 ,5000006 + 02 ,.500000 6+03 ,2500006 + 03 

6 ,5 00000 6 + 0 3 ,.230000E+_03 ^,.7.72 6A.5E+02 .,.5 0 0 0 0 0 6 +.0.3. 

11 .._ . ,9084936 + 02 ,5000006 + 03 ,2500006 + 03 *10.40 046 + 03. 

- 16 . . . ,2500006 + 03 .1171586 + 03 .5000006+03 ,2500006 + 03 

21 , 5000 00 E +0 3 ^,.25_0Q0 06.+i)3 ^»_1 4 50 04.6 + 0 3 ..500000 6 + 0 3 

26 ,15670 2 6 + 03 ^,.500000 6.10 3 , 250000 6 + 03 ,.16 7487 6+03 

31 ,2500006 + 03 , 1774306 + 03. , 5.000006 + 03 .2500006 + 03 

36. ^,.500000 6.+ 03 ,.25^)0006 + 0 3 U 967.7,3 6 + 03 ,.500000 6 + 03. 

41 205 764 6 + 03 ,5000006 + 03 ,.2500006 + 03 ,2141766+03 

46 .2 50000 6+0 3 ,.2 22 2 76 6 + 0 3 , 500000 6. + 03 ^.250000 6 + 0 3 

51 .,5000006 + 03 




NONLINEAR ANALYSIS ITERATION NUMBER 3 


NODE NO. 

NO VALUE 

NO+1 VALUE 

NQf2 VALUE 

NO+3 VALUE 

ND+9 VALUF/ 


1 • 

.250000E+03 

.500000E+02 

. 500000E+03 

25000QF + 03 

. 637703E + 02 

6 

-.500000E + 03 

«250000E+03 

. 775219E+02 

. 500000F+03 

. 250000F+03 

ii 

16 

21 

.913021E + 02 . 
.250000E + 03 
. 500000E + 03 ~ 

.500000E + 03... 
. «119391E+03 
,250000E+03 

.250000E+03. 

. 500000E+03 
. 199R95E + 03 

» 109910E+03 

.250000E+03 

. 500000F+n3 

. ,500000E + 03 

.132569E+03 

. 250000E+03 

26 

. 156993E + 03 

.500000E+03 

.250000F+n3 

. 167999F+03 

» 500000E + 03 

31 

.250000E+03 

.177872E+03 

. 500000F+fl3 

.250000F+03 

.187821F+03 

36 

, 500000E+03 

.250000E+03 

. 197199F ^03 

. 500000F-f03 

•250000E+03 

91 

•206109E+03 

.500000E+03 

. 250000E+03 

» 219987F + 03 

.500000E + 03" 

< 230103E<-03 

96 

.250000E4-03 

».222561E+03 

. 500000E+03 

.250000E+03 

51 

. 500000E+03 





_ 







« ^ _ 

.. . . . 




, 


LARGEST CHANGE . ,,10.2.92O6E+.OI„PE.RCENJ_AT_J^ODE 17. 
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R A_T..U„ R..£ V F C T n R 


„ N 0 N L I N E A R A N AX.Y S I.S_ J1 EBLAJJ.OmJn Um.ER. <l 


NODE NO. NO.__. VAm.e N.0X.l™.m.U.E N.0JL2_y_AXUi Nn+3 VAUJF K.0 + .4__V.ALUEi 


1 • .230000E + 03 ...500000E + 02 . 500000E + 03 .250000E + 03 .637706E + Q2 

.. 6 - .500000Ef03 .250000E.+ 03 ^.J75 26>E+.02 _,.500000E + 03 ,2t)0000E + 03 

11 .913209E + 02 .500000E + 03 .250000E + 03 -104950E + 03 .500000E + 03 

16 . 250000 e + 03 .1193 76E + 03 500000E + 03 .25OQ00E + 03 .132606E + 03 

21 .500000E + 03 ^,Z5Q00 0E+.Q.3 .U 4-4 9 3.Q 6 + 03 ,,500000 E+03. .250000E + 03 

— 26 . _ .156475E + 03 ,500000E+„03 ,.250000E.+.0 3. . 167530E + 03 .500000E + 03 

- 31 .250000E + 03 .177901E + 03 .500000E + 03 .250000E + 03 .187049E + 03 

36 . 50 0000E + 03 ,„Z5.0-0-OOEtO3 ^19.7171E + 03 5000 00E+0 3 2 50000E + 03, 

.. . 41. . .206129E + 03 .500000E + 03 .250000E + 03 .2145116+03 .500000E + 03 

46 .250000E + 03 .222584E+03 .500000E + 03 .250000E + 03 .230124E + 03 

51 .5000 00E+03 


LARGEST CHANGE „,.3.0.7_8_953 E-OI„P„ER.C EN.L_.A_T_N0DE 1 .^ 


C.0 N.V.E R.6 E N.CJ_ACJiL£ALEJl 
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X « G. 0. E.rt.E N 5 I 0 r< A L THE P, HAl-FlU 


li_E A T * FL UX ES 

’TrEHfcTrr'’" fTuid fin condu^tionIfluxes 

HEAT flux’ „QX„_ __ O Y 


1 .1709E + 0A 0* . -,1215E+01 

. ..2 •2I67£ + OA__0.. =.].211>£*.0l 

3 .26A5E+0A 0. -,12L5E+01 

A .3139E»0A 0. -.1215E+01 

. .3663E + 0A _0_, r.»l?rJE + 01 

6 .A2X7E + 0A _.0, =.1215E + 0], 

_ 7 .. .. -A772E+OA__P, r.l215E + 01 

e .S20AE*0 A 0 . !L?_1215E.tOJ 

. , 9 . .i>603E + 0A .0. -.1715E<-01 

10 .59tf3£tO'< .0. 12 15E + 01,_] 

U ,63A!)r + 0^ 0.. .-.l21i»E + 01 

. . - 12 . . .66a8E+0A _0. ^“.121’^£i-0i 

13 _ ,7015E+0A 0. -.1215EF01 

19 ,7326CtO<L_0. T.121&EF01 

15 .7622£+09__0. -.1215E+01 

16 .790AF+09 0. -.1215E*01 



OO OOOjO, o o o o o o 




.0_.V_E._P_ A _L_i L_ I M E LOG 


N 0 D A L . P D I H T __I MP_U .IS 

_F D . E L E y EtU_SJ I .F_F_N.EE:.i_E5_._...,«. , 

F 0 R M T 0 T A L „.S T I EFN E.S.S.,_».._, 

=_ IMPOSE BOUNDARY .CONOIimbiS_, ^.13_ 

E.QUAT ION SOLVING 12_J 

ELEMENT F L UX ES . 1 ? 


.Ti3.LA L_.S ai U T 1,0 N_. LLM E. P.AO 
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SAMPLE PROBLEM 4 


Linear Conduction Analysis of a Fin with Plots 
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INPUT DATA (SAMPLE PROBLEM 4) 


FIN 

PROBLEM CHAPMAN TEXT PAGE 

76. 

EAT NOVEMBER 

14.1 9.75 

1 3 

1 

0 1 

0 1 




1 

1 



0.0 

0.0 

0.0 

2 

0 



1 .0 

0.0 

0.0 

13 

0 



12.0 

0.0 

0.0 

1 

1?_ 






1 

0 

2.75 





1 

1 

2 1 

0 0.19635 

o 

0 

o 

1 .5708 


0.0 

0 

.011111 

70.0 




12 

12 

13 1 

I 0.19635 

0.0 

1 .5708 


0*0 

0 

.011111 

70.0 


■ 



FIN PROBLEM . CHAPMAN TEXT PAGE 76. EAT NOVEMBER 14,1975 
-“t OPT ION 

NNDEST = 13 .NWDISP =1, KPLOT =3. 

£ 

■BPICT 

NOTAT= 1 ,KODE=l ♦ 

S 

apicT 

N0TAT=2.K0DE=1 . 

'I. ' 

SPICT 

NOTAT =l , KDISP =2.DMAGS=0o7. KODE =1, 

•r> 

aPiCT 

NOTAT =2» KDISP =2 . DMAGS=0 « 7 . KODE =1, 
s 

SPICT 

NOTAT=O.KDISP=l . K VERT=3 . KODE= 1 . 
ii» 

aPICT 

KDISP=3.K0DE=0, 

a 



PROGRAM OUTPUT (SAMPLE PROBLEM 4) 


FIM PROflLEH CHAPMAN TEXT PAGE 7 6. 


£AI NQ.V.£tjOER I^t ,I975 


CONTROL INFORHATIHH 


NU MBER OF NODAL PO INTS « 13 


NUMBER OF ELEMENT TYPES - 

1 

' ' 

NUMBER OF TABLES 

0 


ANALYSIS CODE(NANA) - 

1 


EQ.O> DATA CHECK ONLY» 




EQ.1> LINEAR. 


■ 

EQ.Zi NONUNEAR . 

PLOT CODEtNPLOT) « 

1 


EO.O»_NO_PLOTS. generated 



, ^ EQ.l^.. UNOEF0RHED..PJ.0T. 



EQ.2. TEMPERATURE PLOT 

ITERATION PARAMETERS 

- MAXIMUM ITERATIONS 

6 


tolerance 

.lOQOOEtOO 


, 

.lUANK..COrtHON LOCATIONS - 10967 




.NODAL -POINT. INPUT DATA . 




NODE BaUNDARY_CQRD.l.tLDIl_CDILE NODAL pnrN T rno RniWATPC 


NUMBER 


X 

Y 

7 KN 

TEHP.ERATURE 






1 

. 1 

o.noo 

0.000 

0.000 0 

75o*non 

2 

0 

1.000 

0-000 

OfOOO 0 

0.000 

0.000. 

J.3 

0 

IP-.OQO 

0.000 

0.000- 1 


J5E.NE&41Ei)-N0 PAL DATA 


NODE 

BaUNDARY_C0NDlT10K CODE 

NODAL 

filDIHT COORDINATES 


- — ■* - — , . 



X 

Y 

7 KN 

TPMPPRATIIR^= 







1 . _ 

- - 1 

. . 0.000_ 

. .. 0.000 

0.000 

250.000 

0- ooo 

2 

_P 

1.000 

o.noo 

0.000 

3 _ 

.0 . _ 

_ .. __2.000 

0.000 

0 , 000 

. 0.000 

. 0.000 

0.000 
0.000 
0-000 

. ._.A - 

0 

3.000 

0.000 

0-000 

5 

0 

A. 000 

0.000 

0.000 

— 6 „. . 

. .0 

5-000 

n.noo 

0-000 

7 

D - 

.6.000 

0.000 

0.000 

0 

Q 

7.000 

0.000 

0.000 

0 . oon 

5. . 

0 

8.000 

0.000 

0-000 

0 1 000 

10 

0 

9.000 

0.000 

0.000 

0.000 

..11 

0 

10.000 

0.000 

0.000 

0.000 
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— 12 .0 11 *.0.0.0 0,.OCO 0.,J[)0.0 0*000 

13 .0 . „_1Z..000 ._0,000.. 0.000 0.000 



ONE DIMENSIONAL ROD ELEMENT 








NUMBER 

Or ROD ELEMENTS 

= 12 

NUMBER 

OF MATERIALS 

=» 1 


MATERI AL 

CONDUCTIVITY 

CONDUCTIVITY 


TABLE 

K 


1 

0 

.2750E+01 
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CDH0.UC_U.CLM VOLUME 

_ jUS.EA.C2. _ _ 

.Cni<VEC.TJOM. - 

. CDNVEfTIQN DATA 


_ H 

I . 

J 

...HAT 

jvaEA.. 

0 


_A!?£A 

_H1 

TI. 


_ TJ _ 













_ 1 

. 1 

2, 

i ..... 

... ..1969£*00 

0. 

ft. 

.177:F + 01..„ 

.Ullt-01_„ 

. 7000E+02 

.UllE-01 

..70006+02 

2 

2 

3 

. 1 

. .1969E + 00 . 


Qu 

,li71Ci-0I-.. 

.UUE-01 

..7000C»:O2 

.llllErOl. 

-7000E + 02 

3 

3 

9 

1 

,1969e+00 

0. . _ . 

..0 . 

.157U + 01 . 

.UllE-Ol 

•7000E+02 

.llllE-01 

_ .700CE + 0? . . _ 

. S 

- 9 .. 

.5 

« 

.i96't£t00. 

0 . 

0, 

.1571>=4-01 ... 

.niie-oi__ 

.7000E+02. , 

.,11116-01 

•7000F+0? 

> 

5 

6 

1 

.1969E+00 

0. 

.0. 

.lljVlc + Ol 

.llUE-Ol 

.7000E+02 

.llllE-Ol 

. .7000£*02 

6 

6 

7 

1 

.1969£*'00 

0. . 

^0 • _ _ ^ 

.1571E+01 

.llllG-Ol 

.7000£f02 .. 

.11116-01 

. .7000E + 02 .. .. 

. 7 

. 7 

8 

1 

...•1969E + 00_ 

.0. - 


.li71E+01 

.llllE-Ol 

•7000C+32 

..UUE-01. 

70006 + 02 

8 

8 

9 

1 

.1969£f00 

0. 

... ft ^ 

.J571E+01 

.llllE-Ol 

.7000Et02 „ 

.UUE-01 

.70006 + 02 

9 

9 

10 

1 

.1969E+00 

0. 

^ 

.15/lC + Ol 

.iiiir-oi 

.7000E<02 ._ 

UUlE-01 

„ .7000E+02 

. 10 

... 10 . 

11 

_ 1 ._ 

,1969CtO0. 

-0, „ 

0. 

. 157Ut01 

.UllL-Ol _.. 

.7000E*02_. . 

.Ulie-01 

__,7000E+02 

11 

. 11 _ 

12 

.1 

.1969E+00 

.0. 

0 . 

.15/lEt-Ol 

.llllE-Ol 

.7000E402 

.llllc-Ol. 

.7000E+02 

12 

12 

13 

1 - 

. .19C.9E+00 

.0, 

0. . 

. lt> UC + Ol . . 

.llllE-Oi. 

.7000Ef02 

,11116-01 

. ,„.7000E+02 

- 

- 



- 

— 

— 

— 

, , 

— 

. . _ 

— „ 


N) 

CrJ 



$ 0 L U T I □ N 


PARAMETERS 








-TOTAL NUMBER OF EQUATIONS 

ss 

13 

SEMI BANDWIDTH 

Si 

2 

NUMBER OF EQUATIONS IN A BLOCK 

it 

13 

NUMBER OF BLOCKS 


1 
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P_ E. U_E. _E 


V F C T n R 


.NODE N 0 . NO L.U E NO + 1 V.AJJiE _N 04-? V ALUF NO + 3 V At.UF KO + 4 .U E j 


1 .250000E + 03 ,221199E + 03 *19X3l.2E + 03 ♦ 177562E + 03. ♦ 161307E + 03 

- - 6 .148020E+.03. *_X,3.7.2 6 8 E+ Q3 2.8J10.3 E.t.0 3 U2.E.QA5E + 03 1 117079E + 03_ 

11 .. *1136A26 + 03 .11162AE + 03 .n09ij8E + 03 _ , 



0 N E - D : 

[MENSIONA 

L ROD ELEMENTS 


ELEMENT 

CONDUCTION 

SURFACE CONVECTION 


FLUX 

FLUX 



1 

. 15551E+02 

~«28902E+01 

2 

.12898E+02 

~.243GAE+01 

3 

.10664E+02 

-.20A96E^0I 

4 

.87768E+01 

-.1735AE+01 

5 

.71746E+01 

-.1A776E+01 

. 6 

.58055E+01 

-.12679E+01 

7 

.A6251E+01 

10993E+01 

8 

.35950E+01 

-.966A5E+00 

<5 

.268I7E+01 

-.86501 E+00 

10 

.16556E+01 

-.79168E+00 

11 

.10898E+01 

-.7AAC8E+00 

IP 

.35936E+00 

-.72066E+00 




126 





OVERALL 


T I E LOG 


NOD AL . P 0 1 N.T .J N PUJCj. 

^FO_RM_ E.L E NENT S TIF FNESSES..,,.,, ,_U_ 

__FORf1_. TOTAL. STIf1fnESS_ 

IMPOSE 80UN0ARY..C0NDIJ_I_DNS„._._.^. ♦_01_ 

EQUATION. SOLVING . . 1__ 

_ ELEMENT FLUXES, , *D3_ 


XO T A L ..S 0 L U.T 1 0 N _JfX.M E L2A 
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12S 



. 2 3 4 ' 5 6 _7 8 3 ]_Q ll_ 12 13 

(a) Nodes Numbered. 



(b) Elements Numbered. 


5 6 7 8 3 10 y 12 13 

(c) Exploded with Nodes Numbered. 

_p _5 _g _Q _.o _4^3_ 

(d) Exploded with' Elements Numbered. 

Figure 16. Plotter output for sample problem 4. 
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